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^ ■ Abstract 

The SU(Nf)L® SU(Nf)n chiral symmetry of QCD is of central importance for the nonpertur- 
bative low-energy dynamics of light quarks and gluons. Lattice field theory provides a theoretical 
7— I ' framework in which these dynamics can be studied from first principles. The implementation of 

chiral symmetry on the lattice is a nontrivial issue. In particular, local lattice fermion actions 
| with the chiral symmetry of the continuum theory suffer from the fermion doubling problem. The 

Ginsparg- Wilson relation implies Liischer's lattice variant of chiral symmetry which agrees with 
the usual one in the continuum limit. Local lattice fermion actions that obey the Ginsparg- Wilson 
relation have an exact chiral symmetry, the correct axial anomaly, they obey a lattice version of 
the Atiyah-Singer index theorem, and still they do not suffer from the notorious doubling problem. 
The Ginsparg- Wilson relation is satisfied exactly by Neuberger's overlap fermions which are a limit 
of Kaplan's domain wall fermions, as well as by Hasenfratz and Niedermayer's classically perfect 
lattice fermion actions. When chiral symmetry is nonlinear ly realized in effective field theories on 
the lattice, the doubling problem again does not arise. This review provides an introduction to 
chiral symmetry on the lattice with an emphasis on the basic theoretical framework. 
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1 Introduction 



Physical phenomena arise over a vast range of energy scales. Attempts to unify gravity with the other 
fundamental forces suggest that the Planck scale 

M P = ^=^ 10 19 GeV, (1) 
VG 



which is constructed from Newton's constant G (and from % and c which we have put to 1) is the highest 
energy scale relevant to particle physics. On the other hand, ordinary matter receives almost all of its 
mass from protons and neutrons which have a mass M « 1 GeV. Can we understand why nucleons 
exist so far below the Planck scale? This is a typical hierarchy problem of which there are several in 
physics (also including the notorious cosmological constant problem). As Wilczek has pointed out, the 
large hierarchy between the Planck scale Mp and the nucleon mass M is responsible for the feebleness 
of gravity PQ. To understand this, let us compare the strengths of the gravitational attraction and the 
electrostatic repulsion between two protons at some large distance R. The force of gravity is given by 

M 2 

F 9 = G w , (2) 



while the electrostatic force is 



e 2 



(3) 



where e is the proton's electric charge. The ratio of the two forces is thus 



§ = G% « !37^ „ 10- (4) 



e 



p 



Hence, if we can understand why M <C Mp, we can understand why gravity is a very weak force. 

As Wilczek has explained in pQ, the nucleon mass M is much smaller than the Planck scale Mp 
partly due to the asymptotic freedom of QCD. At the classical level, QCD with massless quarks has 
no dimensionful parameter at all. When the theory is quantized, a scale enters through the mecha- 
nism of dimensional transmutation. A non-zero nucleon mass arises even in massless QCD due to the 
confinement of colored quarks and gluons inside color-neutral hadrons. Thus, the nucleon mass M is 
a nonperturbatively generated scale which cannot be understood using perturbation theory. The con- 
tinuum theory, i.e. dimensional regularization and renormalization applied to the QCD Lagrangian, is 
not even defined at a nonperturbative level. The only nonperturbative definition of QCD is provided 
by lattice field theory in which continuous space-time is replaced by a hypercubic lattice [21 EJ- In 
particular, it should be pointed out that lattice QCD is not an approximation to any pre-existing non- 
perturbatively well-defined theory in the continuum. Of course, as in any other quantum field theory, 
one must ultimately remove the cut-off. On the lattice, the shortest physical distance is the lattice 
spacing a which defines an ultraviolet momentum cut-off 1/a. Removing the cut-off thus means taking 
the continuum limit a — > 0. The masses of hadrons M — l/£ are the inverse of a correlation length £. 
Taking the continuum limit means that the physical mass M must be much smaller than the cut-off, 
i.e. 

M < - i > a. (5) 
a 

Hence, in the continuum limit the physical correlation length £ goes to infinity in units of the lattice 
spacing. In the language of classical statistical mechanics this corresponds to a second order phase 
transition. 

Most of the time, lattice QCD is used as a very powerful tool for numerical calculations of hadronic 
properties. However, the lattice can do more than that. To illustrate this, we will now use lattice 
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QCD to explain why nucleons can exist naturally far below the Planck scale. Of course, it is well- 
known that QCD is not valid up to the Planck scale. In particular, it is embedded in the standard 
model which itself is an effective theory limited to energies below about 1 TeV. However, unlike the 
full standard model, thanks to asymptotic freedom QCD alone makes sense at arbitrarily high energy 
scales. Whatever replaces QCD and the standard model at ultra-short distances — be it string theory 
or some tiny wheels turning around at the Planck scale — Nature must have found a concrete way to 
regularize the QCD physics at ultra-short distances. Due to renormalizability and universality, only 
the symmetries but not the details of this regularization should matter at low energies. For simplicity, 
we will use lattice QCD (and not, for example, string theory) as an admittedly oversimplified model of 
Nature at ultra-short distances. In other words, in this context we like to identify the lattice cut-off 
l/a with the Planck scale Mp. 

Using lattice QCD, how can we then understand why the nucleon mass M is far below Mp = l/a? 
As Wilczek pointed out, one key ingredient is asymptotic freedom. Thanks to asymptotic freedom, 
without any fine-tuning of the bare gauge coupling a non-Abelian lattice Yang-Mills theory produces a 
correlation length £ that is larger than the lattice spacing a by a factor exponentially large in the inverse 
coupling. In particular, choosing a bare coupling that is not unnaturally small, one can easily generate 
a hierarchy of scales like Mp/M = £/a ~ 10 19 . Interestingly, the situation is not as simple when one 
proceeds from a pure gluon Yang-Mills theory to full lattice QCD including quarks. In particular, unlike 
continuum QCD, lattice QCD does not naturally have a chiral symmetry that can protect the quark 
masses from running up to the lattice "Planck scale" l/a. Indeed, for about two decades lattice field 
theorists have suffered from a hierarchy problem in the fermion sector. This problem first arose when 
Wilson removed the unwanted doubler fermions by breaking chiral symmetry explicitly [Sj. Recovering 
chiral symmetry in the continuum limit then requires a delicate fine-tuning of the bare fermion mass. 
In particular, if at ultra-short distances Nature would be a lattice gauge theory with Wilson fermions, 
without unnatural fine-tuning quarks would have masses at the Planck scale and the lightest particles 
would be glueballs. In that case it would be very puzzling why ordinary matter consists not just of 
gluons, but also of light quarks. If one works in continuum QCD one often takes chiral symmetry for 
granted, and one may view this hierarchy puzzle just as a problem of the lattice formulation. However, 
one should not forget that continuum QCD is not even defined beyond perturbation theory. In addition, 
subtleties of the definition of 75 in the framework of dimensional regularization affect even the continuum 
theory, and are just another aspect of the same deep problem of chiral symmetry that is manifest on the 
lattice. Indeed, there is a severe hierarchy problem for nonperturbative fermion dynamics that Nature 
must have solved somehow because it presents us with nucleons that exist far below the Planck scale. 

Remarkably, the long-standing hierarchy problem of lattice fermions — and hence of the nonper- 
turbative regularization of chiral symmetry — has recently been solved very elegantly. Using previous 
research of Callan and Harvey jl], Kaplan [3] realized that massless four- dimensional lattice fermions 
arise naturally, i.e. without fine-tuning, as zero-modes localized on a domain wall embedded in a five- 
dimensional space-time. In particular, left- and right-handed fermions can be localized on a domain 
wall and an anti-wall. When the wall and the anti-wall are separated by a sufficiently large distance, 
the left- and right-handed modes cannot mix, simply because they are spatially separated. As a result, 
a Dirac fermion arises which is protected from picking up a large mass and which is thus naturally light. 
Remarkably, in contrast to four dimensions, a Wilson term in a five-dimensional lattice theory removes 
the doubler fermions without breaking the chiral symmetry of the light four- dimensional domain wall 
fermions. 

When Kaplan proposed his idea of regulating chiral fermions using domain walls, Narayanan and 
Neuberger were developing independently another approach to regulating chiral fermions using an infi- 
nite number of flavors [6 . Based on this approach they developed what is now referred to as overlap 
lattice fermions [HE]. Since the flavor-space can be viewed as an extra dimension, the overlap approach 
is closely related to the domain wall approach. When one separates the wall and the anti-wall by an 
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infinite distance, domain wall fermions turn into overlap fermions. Overlap fermions have the advan- 
tage that they have an exact chiral symmetry, while the chiral symmetry of domain wall fermions is 
only approximate for a finite wall-anti-wall separation. Both overlap and domain wall fermions yield 
naturally light quarks, and both are naturally related to the physics of a higher-dimensional space-time. 

Hasenfratz and Niedermayer have investigated nonperturbative renormalization group blocking 
transformations on the lattice [U] . The fixed points of such transformations correspond to lattice actions 
which are completely free of cut-off effects — so-called perfect actions. Perfect actions for pure gauge 
theory, as well as for free Wilson and staggered fermions were investigated in piHITHIT21 ll3j. Classically 
perfect actions for full QCD have been constructed in [TU[T3]. In the process of these investigations 
Hasenfratz (see [THj) rediscovered an old paper by Ginsparg and Wilson [T7]. He also realized that what 
is now called the Ginsparg- Wilson relation is the key to understanding chiral symmetry on the lattice. 
The Ginsparg- Wilson relation represents a general requirement on a lattice action which guarantees 
that it has good chiral properties. When Ginsparg and Wilson discovered this relation, it seemed im- 
possible to explicitly construct lattice actions that obey it. By now it has been shown that classically 
perfect lattice actions can be approximated well enough, so that the Ginsparg- Wilson action is satis- 
fied with high accuracy. From the point of view of practical lattice QCD calculations this represents 
very important progress. However, the explicit construction of perfect actions is a delicate problem 
that can be considered a very elaborate form of fine-tuning. Hence, it seems unnatural that Nature 
has chosen anything like a perfect action to regularize the strong interactions at ultra-short distances. 
Unlike perfect fermions, overlap fermions can describe massless quarks in QCD without fine-tuning. 
By integrating out the extra dimension, Neuberger has constructed lattice Dirac operators for massless 
quarks analytically fHl EHj and these Dirac operators do indeed satisfy the Ginsparg- Wilson relation 
exactly [201 • Remarkably, both overlap as well as domain wall fermions, which naturally have a chiral 
symmetry without fine-tuning, are related to the physics in a higher- dimensional space-time. Hence, 
the existence of light four-dimensional fermions may be a concrete hint to the physical reality of extra 
dimensions. 

The full strength of the Ginsparg- Wilson relation was realized by Luscher who discovered that it 
suggests a natural definition of lattice chiral symmetry which reduces to the usual one in the continuum 
limit [21]. Based on this insight, Luscher achieved a spectacular breakthrough: the nonperturbative 
construction of lattice chiral gauge theories [221- Hence, not only QCD in which chiral symmetry is 
global, but also the standard model with its local chiral symmetry now stands on a solid nonperturbative 
basis. Even continuum perturbation theory can benefit from these developments. In particular, the 
ambiguities in the definition of 75 that arise in multi-loop calculations using dimensional regularization 
can be eliminated when one uses the lattice regularization. Still, there is a very long way to go from 
Liischer's theoretical construction to practical numerical calculations in chiral gauge theories like the 
standard model. 

The situation is a lot simpler, but still highly nontrivial, in applications of Ginsparg- Wilson fermions 
to simulations of QCD. Compared to the standard Wilson or staggered lattice fermions, which are 
already very difficult to treat fully dynamically, domain wall, overlap, or perfect fermions demand even 
much larger computing power. Hence, at present they are often used in the quenched approximation 
in which the fermion determinant is ignored. If one does not want to wait a long time for even bigger 
computers, it will require an algorithmic breakthrough to bring the theoretical developments of lattice 
chiral symmetry to fruition in fully dynamical simulations of lattice QCD. Promising steps in this 
direction are already being taken |23] [21] . 

If one imagines that Nature has used something like domain wall fermions to regularize the strong 
interactions, it is natural that nucleons (and not just glueballs) exist far below the Planck scale. How- 
ever, it remains mysterious where the quark masses themselves come from. In the standard model the 
quark masses arise from Yukawa couplings to the Higgs field, but the values of these couplings are free 
parameters. Still, within the standard model the traditional gauge hierarchy problem arises: why is the 
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electroweak scale so small compared to the Planck scale? Chiral symmetry can protect fermion masses 
from running to the ultimate high-energy cut-off, but it cannot protect scalars. A potential solution of 
the gauge hierarchy problem is provided by supersymmetry. Supersymmetry relates scalars to fermions 
and thus allows chiral symmetry to indirectly generate naturally light scalars as well. At a nonper- 
turbative level, supersymmetry is as undefined as chiral symmetry was before the recent developments 
on the lattice. In the worst case, supersymmetry may just be a perturbative illusion which does not 
arise naturally at a nonperturbative level. Unfortunately, unlike for chiral symmetry, Nature has not 
yet provided us with experimental evidence for supersymmetry (except as an accidental symmetry in 
heavy nuclei). Hence, one cannot be sure that it is indeed possible to construct naturally light scalars 
at a nonperturbative level. Perhaps the many beautiful results obtained within supersymmetric con- 
tinuum theories should make us optimistic that these theories actually exist at a rigorous level beyond 
perturbation theory. Again, Kaplan and his collaborators have taken very interesting steps towards 
constructing supersymmetric theories on the lattice [25J. It remains to be seen if these developments 
will lead to a repetition of the Ginsparg- Wilson revolution of lattice chiral symmetry. 

This review is an introduction to chiral symmetry on the lattice. We assume that the reader has 
a background in continuum field theory but not necessarily in lattice field theory. It should be noted 
that there are very good textbooks that cover the basics of lattice field theory much more completely 
than it can be done here [211 123 EH] ■ Furthermore, there are excellent reviews of the more advanced 
and specialized aspects of lattice field theory and chiral symmetry |29} l30\ l3Tj 132] . This review aims at 
bridging some of the gaps between the basic and the more advanced texts and introducing the reader to 
some of the latest developments in the field. In section 2 we summarize symmetry properties of QCD in 
the continuum. We concentrate on chiral symmetry, but also discuss scale invariance, the axial anomaly, 
and the Atiyah-Singer index theorem. Section 3 illustrates basic properties of lattice chiral symmetry 
using free fermions. In particular, we discuss the fermion doubling problem, the Nielsen-Ninomiya 
theorem, and Wilson, staggered, as well as perfect fermions. In section 4, gluons are added and full 
lattice QCD is considered. In particular, the Ginsparg- Wilson relation is discussed. Section 5 concerns 
special features of Ginsparg- Wilson fermions, including the axial anomaly and the index theorem on 
the lattice, the Witten-Veneziano mass formula, the renormalization of operators, as well as numerical 
simulations of Ginsparg- Wilson fermions. Section 6 discusses effective theories for pions, nucleons, and 
constituent quarks, both in the continuum and on the lattice. In particular, the nonlinear realization 
of chiral symmetry implies another way of solving the fermion doubling problem. Lattice simulations 
of constituent quarks may eventually shed some light on the success of the nonrelativistic quark model. 
Finally, section 7 contains our conclusions. 



2 Symmetries of the Strong Interactions 

In this section we review some aspects of symmetries in the continuum formulation of QCD with an 
emphasis on chiral symmetry. 

2.1 SU(N C ) Yang-Mills Theory in the Continuum 
Let us consider an anti-Hermitean non-Abelian SU(N C ) gauge field 

A,(x) = igA°(x)T a , (6) 

which (for N c = 3) describes the gluons of QCD. Here g is the gauge coupling, A a Jx) (with a G 
{1,2,..., N% — 1}) is the real- valued non-Abelian vector potential at the Euclidean space-time point x, 
and the T a (which obey Tr(T a T b ) = \5 a b) are the Hermitean generators of the SU(N C ) algebra. The 
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algebra-valued field strength takes the form 

F^(x) = d»A v {x) - d u A^(x) + [A^x), A„(x)], (7) 
and the corresponding Euclidean Yang-Mills action is given by 

Sym[A] = -Jd 4 x -lTr(F^). (8) 

The action is invariant under group- valued gauge transformations Q(x) G SU(N C ), 

A' ll (x) = n(x)(A ll (x) + d lt )n(x)i, (9) 
under which the field strength transforms as 

f^(x) = n(x)F^(x)n( x y. (io) 

The quantum theory is defined by a functional integral over all gluon fields 

Z = JVA exp(-S YM [A}), (11) 

which is a formal expression before it is properly regularized. In perturbation theory this is possible 
using standard dimensional regularization techniques. Through the regularization, a scale is introduced 
into the quantum theory which explicitly breaks the scale invariance of the classical Yang-Mills theory. 
This anomaly in the scale invariance is responsible for the phenomenon of dimensional transmutation: 
in the quantum theory the dimensionless coupling constant g of the classical theory is traded for a 
dimensionful scale. In the modified minimal subtraction renormalization scheme this scale is 
which is defined in the framework of perturbation theory. We will soon define the theory beyond 
perturbation theory by regularizing it on a space-time lattice. In a nonperturbative context, a natural 
scale is the dynamically generated mass gap M — the energy of the lowest state above the vacuum. 
In a Yang-Mills theory this state is the lightest glueball. The SU (N c ) Yang-Mills theory is a quantum 
theory without any free parameter. For example, the dimensionless ratio M/Ajjg is a pure number 
predicted by the theory. The relation of M or Aj^ to units like GeV, on the other hand, is, of course, 
not predicted by the theory. Such man-made mass units are related to the kilogram, defined by the 
arbitrary amount of platinum-iridium alloy deposited near Paris a long time ago. 
Another quantity of physical interest is the topological charge 

Q[A] = J d 4 x e flvper Tr(F fU/ F pa ) G U 3 [SU(N C )} = Z, (12) 

which takes integer values in the third homotopy group of the gauge group. The topological charge 
gives rise to an additional parameter, the vacuum angle 9, in the Yang-Mills functional integral 

Z{9) = JVA exp(-S Y M[A]+i9Q[A]). (13) 

For 6 7^ or ir the #-term explicitly breaks parity as well as CP. The bound \9\ < 1CT 9 derived from the 
measurement of the electric dipole moment of the neutron suggests that 9 = in Nature. This result is 
puzzling because in the Standard Model CP is already explicitly broken by the complex phase of the 
Cabbibo-Kobayashi-Maskawa matrix. The puzzle to understand why 6 = is known as the strong CP 
problem. 
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2.2 QCD with N f Quark Flavors 

In the next step we add Nf massless quarks to the pure gluon theory. Quarks and anti-quarks are 
described by anti-commuting Dirac spinor fields ifj(x) and ip(x). In Euclidean space-time these two fields 
represent independent Grassmann degrees of freedom. Under a non-Abelian gauge transformation the 
quark and anti-quark fields transform in the fundamental representations {iV c } and {N c }, respectively, 
i.e. 

i>(x)' = n(x)ip(x), $(x)' = ^(x)n(x)l (14) 

The fermionic part of the Euclidean action of massless QCD takes the form 



S, 



#,^M] = J d^^d^ + A^, (15) 



which is gauge invariant by construction. The Euclidean Dirac matrices are Hermitean and obey the 
anti-commutation relations 

{l^lu} = 2<V, {7^75} = 0, 75 = 7i727s74- (16) 

We now decompose the quark fields into left- and right-handed components 

ip L (x) = P L il>(x), ip R (x) = P R ip(x), i/j(x) = ip L (x) +i/j r (x), 

tf L (x) = Hx)P R , tf R (x) = tf(x)P L , ^(x) = i> L (x) + tf R (x). (17) 



The chiral projectors are given by 



ft = i±* ft = i^i (18) 



Inserting the decomposed spinors into the fermionic part of the action one obtains 

S F ®, ip,A] = J d 4 x p L7/i (^ + + ^7m(^ + a Mr] > ( 19 ) 

i.e. the action decouples into two contributions from left- and right-handed quarks. 

As a result, the action of massless QCD is invariant against U (Nf) L ®U(Nf) R chiral transformations 

tf L (x) = L M*), = ^l(x)L + , L E U(N f ) L , 

i/>' R (x) = R M%), ^(x) = ^ R (x)R + , R e U(N f ) R . (20) 

Due to an anomaly in the axial U{1)a symmetry, the symmetry of the quantum theory is reduced 
to SU(Nf) L (g> SU(Nf) R <g> U(1)b where the U(1)b = U(\)l=r symmetry represents baryon number 
conservation. 

Chiral symmetry is only approximate in Nature, because the quark mass terms couple left- and 
right-handed fermions. The mass terms in the QCD action take the form 

S M ®,il>] = J d A x \4> R Mi> L + ^ L MH R ] , (21) 

which is again gauge invariant but no longer chirally invariant. The quark mass matrix takes the form 

M = diag(m M , m d , m s , m Nf ). (22) 

If all quark masses are equal, i.e. if M. = ml, the mass term is invariant only against simultaneous 
transformations L = R. Hence, chiral symmetry is then explicitly broken down to 

SU(N f ) L=R <g> U(1) L=R = SU{N f ) F <g> U(1) B , (23) 
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which corresponds to the flavor and baryon number symmetry. In Nature the quark masses are all 
different, and the symmetry is, in fact, explicitly broken down to 

J] U(l) f = U(l) u ® U(l) d ® U(l) 3 ® ... ® U(l) Nf . (24) 
/=i 

The physical up and down quark masses are a lot smaller than while the strange quark mass is 
of the order of A^. Consequently, SU(2)l ® SU(2)r is a very good approximate global symmetry, 
while SU(3)l <8> SU(3)r is broken more strongly It should be noted that the actual values of the 
quark masses are reasonably well-known from comparison with experiment, but are at present not at all 
understood theoretically. In particular, we don't know why there are three light quark flavors. Before 
one understands the relevant physics beyond the Standard Model, the origin of the chiral symmetry of 
QCD remains mysterious and the symmetry itself seems accidental. 
The total action of QCD is simply given by 

S QCD $i V, A] = S YM [A] + S F $, V, A] + S M ®, V], (25) 
and the corresponding QCD functional integral is 

Z = J ThJtihpVA exp(-S QCD $, <M]). (26) 

Again, this is a formal mathematical expression before it is properly regularized. In the continuum this 
can be done only perturbatively. We will soon discuss the lattice regularization which defines QCD 
beyond perturbation theory. 

2.3 The Axial Anomaly and the Atiyah-Singer Index Theorem 

The U(1)a symmetry of the classical action of massless QCD is explicitly broken by quantum effects. 
As a consequence of this anomaly the flavor- singlet axial current 

jl(x) =^(x) 7m 7 5 ^(x), (27) 

which is conserved at the classical level, has a non-zero divergence 

d^x) = -^ W Tr[F^(x)F pa (x)], (28) 

due to instantons (and other topological charge carriers) in the quantum theory. In particular, the 
variation of the axial charge Q 5 (t) = J d 3 x jo(x,t) is given by 

Q 5 (t = oo)-Q 5 (t = -oo) = N f Q[A], (29) 

where Q[A] is the topological charge of eq. (jl2j) . Only in the N c —>■ oo limit the anomaly vanishes and 
the chiral symmetry of massless QCD is enhanced to the full U(Nf)i ® U(Nf)n group. 

The axial anomaly is deeply connected with the Atiyah-Singer index theorem, which relates the 
zero- modes of the massless (iVy-flavor) Dirac operator D[A] = 7 M (i9 M + A^) to the topological charge 
Q[A]. The eigenvalues of the Dirac operator are purely imaginary and come in complex conjugate 
pairs. Only the zero eigenvalues are not paired. Since the Dirac operator anti-commutes with 75, the 
eigenvectors of the zero-modes (which obey DL4]^ = 0) have a definite handedness, i.e. 75^ = 
The index theorem states that 

index(D[A]) = n_ - n + = N f Q[A], (30) 
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i.e. the index of the operator -D[A], which is defined as the difference between the number of left- and 
right-handed zero-modes, is given by the topological charge. 

As a consequence of the index theorem, topologically nontrivial gluon field configurations (with 
Q[A] 7^ 0) necessarily induce zero- modes in the Dirac operator and thus lead to a vanishing fermion 
determinant det/}[7L] = 0. As a function of the vacuum angle 9, the functional integral of massless 
QCD takes the form 

Z{9) = J V^VtPVA exp(-S QCD $, V>, A] + i9Q[A\) 
= jvA detD[A] exp(-S YM [A} + i9Q[A}) 

= JvA detD[A] exp(-£y M L4]) = Z(0). (31) 

Since detDL4] = when Q[A] ^ 0, there are no vacuum effects in massless QCD. This would "solve" 
the strong CP problem (why is 9 = ?) if, for example, the up quark would be massless. Of course, 
this would leave us with the "up quark problem" : why should m u = ? In any case, m u = seems not 
to be realized in Nature and the strong CP problem remains puzzling. 

2.4 Spontaneous Chiral Symmetry Breaking 

Due to the approximate chiral symmetry of QCD one would expect corresponding near degenera- 
cies in the spectrum of strongly interacting particles. Indeed, hadrons can be classified as isospin 
multiplets. The isospin transformations act on left- and right-handed fermions simultaneously, i.e. 
£77(2)/ = SU(2) L=R . The SU(3)p = SU(3)l =r flavor symmetry is more approximate but is still 
clearly visible in the spectrum. The full SU(Nf)i ® SU(Nf) R <2)U(1)b chiral symmetry, on the other 
hand, is not manifest in the spectrum at all. In particular, one does not observe mass-degenerate parity 
doublets of hadrons, as one should if chiral symmetry was manifest in the spectrum. Furthermore, one 
observes very light pseudo-scalar particles — the pions n + , n°, and ir~ - as well as somewhat heavier 
pseudo-scalars — the four kaons K + , K°, K°, K~ and the 77- meson. 

/,From the experimental evidence one concludes that chiral symmetry must be spontaneously broken. 
Indeed, when a continuous global symmetry breaks spontaneously, massless Goldstone bosons appear 
in the spectrum. According to Goldstone's theorem, the number of massless bosons is given by the 
difference of the number of generators of the full symmetry group G and the subgroup H that remains 
unbroken. In massless QCD the full chiral symmetry group is 

G = SU(N f ) L ®SU(N f ) R ®U(l) B , (32) 

while the unbroken subgroup is the flavor symmetry 

H = SU(N f ) L=R ®U(l) B . (33) 

Hence, in this case one expects Nj — 1 massless Goldstone bosons. For Nf = 2 these are the three 
pions, while for Nf = 3 there are eight Goldstone bosons — the pions, the kaons, and the r/-meson. 
In Nature these particles are not exactly massless, because chiral symmetry is explicitly broken by the 
quark masses. The masses of the up and down quarks are much smaller than the QCD scale A-p^ which 
leads to the very small pion mass. The mass of the strange quark, on the other hand, is of the order of 
A^r, thus leading to larger masses of the kaons and the ?7-meson. Still, their masses are small enough 
to identify these particles as pseudo-Goldstone bosons. 

Chiral symmetry breaking has not yet been derived analytically from the QCD Lagrangian. In par- 
ticular, spontaneous chiral symmetry breaking is a nonperturbative phenomenon whose understanding 



9 



requires a formulation of QCD beyond perturbation theory. Such a formulation is provided by lattice 
field theory which will be discussed below. In the lattice formulation, in the strong coupling limit it 
is possible to show rigorously that chiral symmetry is indeed spontaneously broken |B*3*1 . Further, 
numerical simulations in lattice QCD confirm that chiral symmetry is spontaneously broken even at 
weaker couplings. For example, one detects spontaneous chiral symmetry breaking by investigating the 
chiral order parameter 

®i/>) = ((#(x)V(a;)|0) = (0® R (x)Mz) +^ L {x)i/) R (x)\0). (34) 

The order parameter is invariant against simultaneous transformations R = L, but not against general 
chiral rotations. If chiral symmetry would be intact the chiral condensate would vanish. When the 
symmetry is spontaneously broken, on the other hand, (ipip) is non-zero. 



3 Free Lattice Fermions 

In this section we begin to formulate QCD on a space-time lattice which serves as an ultraviolet regulator. 
We replace Euclidean space-time by a hypercubic lattice of points x with lattice spacing a. The lattice 
provides an ultraviolet momentum cut-off 1/a. The continuum limit is reached when a — > 0. One has 
a lot of freedom in writing down a lattice regularized theory. In order to ensure that one reaches the 
desired theory in the continuum limit one must pay attention to the relevant symmetries. The most 
important symmetry of QCD is the SU(N C ) gauge invariance. It is an important strength of the lattice 
regularization that it manifestly respects gauge invariance. The fact that space-time symmetries are 
explicitly broken down to discrete translations and the hypercubic rotation group of the lattice is not a 
severe problem. In particular, in QCD the hypercubic symmetry is powerful enough to ensure that the 
full Poincare symmetry of the continuum is automatically recovered as a — > 0. Discrete symmetries like 
parity and charge conjugation are also easy to maintain on the lattice. This review concentrates on the 
question of how to realize chiral symmetry on the lattice. In this section we consider lattice theories of 
free quarks only. Gluon fields will be added in the next section. 



3.1 The Naive Lattice Fermion Action and the Doubling Problem 

In the continuum the Euclidean action of a free Dirac fermion in d space-time dimensions is given by 

S$, V>] = J d d x ^( 7 A + m)V, (35) 

and the functional integral takes the form 

Z = Jv^Vip exp(-5f0,V]). (36) 

On the lattice the continuum fermion field ip(x),ip(x) is replaced by Grassmann variables ty x , ^> x which 
live on the lattice points x. The continuum derivative can be discretized by a finite difference, such that 

tf] = a d VDV = a d J2 ^-(^x7m**+/> ~ **+A7„^) + E (37) 

Here fi is a vector of length a in the //-direction. In the continuum limit a — > the lattice sum a d J2 X 
becomes the continuum integral / d d x over space-time. The corresponding lattice Dirac operator which 
is a matrix in the Dirac- and space-time indices takes the form 
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The lattice functional integral can be written as 

Z = J V^V^ exp(-a d ^D^) =11 J <^ x d^ x exp(-5[*, *]). (39) 

In particular, the fermionic Grassmann integration measure is completely regularized explicitly. 

The momentum space of the lattice theory is a <i-dimensional Brillouin zone B = [—n/a, ix/a\ d with 
periodic boundary conditions. Going to momentum space, the naive fermion action from above gives 
rise to the lattice fermion propagator 

(^(-p)tf(p)) = [^7/,-sinGvi) +m}-\ (40) 
n a 

By performing a Fourier transform in the Euclidean energy pd one obtains the fermion 2-point function 
0)*(p, = r ° (*(-p)*(p)> exp(^x d ) ~ exp(-£(p> d ). (41) 

Z7T J-Tr/a 

At large Euclidean time separation the 2-point function decays exponentially with the energy E(p) 
of a fermion with spatial momentum p. For the naive fermion action the lattice dispersion relation takes 
the form 

smh 2 (E(p)a) = ^ sin 2 (p l a) + (ma) 2 . (42) 

i 

The continuum dispersion relation E(p) 2 = p 2 + m 2 is indeed recovered in the continuum limit a — > 0. 
However, besides p = there are other momenta p for which E(p) becomes small. These are located at 
the corners of the Brillouin zone where the components of the momentum vector take the values Pi = 
or 7i /a, such that sin(pja) = 0. As a consequence, the lattice dispersion relation leads to additional 
states in the spectrum which are absent in the continuum theory and which do not disappear in the 
continuum limit. Hence, the naive lattice fermion action does not lead to the correct continuum theory. 
The extra states appearing in the lattice dispersion relation show up as additional physical particles 
- the so-called doubler fermions. Fermion doubling is a manifestation of a deep fundamental problem 
of lattice regularized fermionic theories with a chiral symmetry. The fermion doubling problem leads 
to a multiplication of fermion species. The lattice fermion propagator of eq. (}4T)|) has 2 d poles instead 
of just one as in the continuum. The origin of the doubling problem is deeply connected with chiral 
symmetry and can be traced back to the axial anomaly. The doubler fermions pose a severe problem 
in lattice field theory. Without removing them we cannot describe Nature's QCD (which has 3 and not 
2 d = 2 4 = 16 light quark flavors). 

3.2 The Nielsen-Ninomiya Theorem 

Before we try to eliminate the doubler fermions let us prove a general theorem due to Nielsen and 
Ninomiya [33] : a chirally invariant free fermion lattice action, which is local, translation invariant, and 
real necessarily has fermion doubling. The theorem is based on topology. It holds because the lattice 
momentum space (the Brillouin zone B) is a torus. A general chirally symmetric and translationally 
invariant lattice action for free fermions takes the form 

=a i Y^xlM^-V)%- ( 43 ) 

x,y 

The function p M (x — y) determines the strength of the coupling between the fermion field values ty x 
and ty y at two points x and y which may be separated by an arbitrarily large distance. Locality of 
the lattice action does not mean that the points x and y must be nearest neighbors. It only means 
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that p^{x — y) decays exponentially at large separations x — y. Going to momentum space, locality 
implies that in Fourier space p^ip) is a regular function (without poles) over the Brillouin zone. The 
corresponding lattice fermion propagator takes the form 



Reality and translation invariance of the lattice action imply that p^ip) is a real- valued periodic function 
over the Brillouin zone. 

Poles of the propagator — and hence physical or doubler fermions — correspond to zeros of p^Qo), 
i.e. to points p with = for all p. The Nielsen- Ninomiya theorem states that a regular, real- valued, 
and periodic function p,j,(p) necessarily vanishes at more than just one point. It is trivial to prove this 
for d — 1. In that case, there is a single regular periodic function Pi(p) which should at least have one 
zero in order to describe the physical fermion pole. The function is positive on one side of the zero and 
negative on the other side. Hence, it must go through zero again in order to satisfy periodicity, thus 
leading to a doubler fermion pole in the lattice propagator. In higher dimensions the proof is analogous. 
For example, for d — 2 there are two functions pi(p) and Pzip)- The zeros of p\{p) lie on a closed curve 
in the two-dimensional Brillouin zone. This curve may be closed via the periodic boundary conditions. 
The zeros of p2(p) he on another closed curve that intersects the first one in the pole position of the 
physical fermion. Due to the periodic boundary conditions of the Brillouin zone, the two curves must 
necessarily also intersect somewhere else. The curves cannot just touch each other because this would 
lead to an incorrect dispersion relation for the physical fermion. In d dimensions the zeros of p^ip) (with 
p = 1,2, ...,d) lie on d closed (d — l)-dimensional surfaces. Again, those cannot intersect in just one 
point. If they intersect once they necessarily intersect also somewhere else. This proves lattice fermion 
doubling for a chirally symmetric, translation invariant, real- valued lattice action. It should be noted 
that the theorem does not specify the number of doubler fermions. It is indeed possible to reduce the 
number of doublers from 2 d — 1 to 1, but it is impossible to eliminate the doubler fermions completely. 

One may try to evade the theorem by violating one of its basic assumptions. Giving up translation 
invariance or the reality of the action has not led to acceptable solutions of the fermion doubling 
problem. Giving up locality is probably the last thing one should do in field theory. For example, the 
early idea of SLAC fermions [SH] turned out to be unacceptable for this reason. 

3.3 Wilson Fermions 

In his work on lattice gauge theory Wilson removed the doubler fermions in a direct and radical way 
by breaking chiral symmetry explicitly 3J. Then the Nielsen- Ninomiya theorem is evaded because the 
propagator contains additional terms without 7 M . The so-called Wilson term gives the fermion doublers 
a mass of the order of the cut-off while the physical fermion remains massless. Hence, in the continuum 
limit chiral symmetry is recovered in the physical sector. Wilson's modification of the naive fermion 
action takes the form of a discretized second derivative 




(44) 




(45) 



Then the lattice propagator takes the form 



sin (— )] • 



(46) 
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The Wilson term acts as a momentum-dependent mass term. For small momenta it vanishes quadrat- 
ically, and hence it does not affect the dispersion of the physical fermion, at least in the continuum 
limit. For the doubler fermions, on the other hand, the Wilson term is non-zero, and gives them a mass 
of the order of the cut-off 1/a. In the continuum limit the doubler fermions are hence eliminated from 
the spectrum of the theory. Unfortunately, in lattice QCD this leads to a variety of complications. In 
particular, recovering chiral symmetry in the continuum limit requires unnatural fine-tuning of the bare 
fermion mass. For more details we refer the reader to [27| I28j. 

3.4 Staggered Fermions 

Staggered fermions result from naive doubled lattice fermions by so-called spin diagonalization [37J EH1 
13*9*] . In four space-time dimensions the size of the Dirac matrices is 4 x 4. By spin diagonalization one 
can reduce the fermion multiplication factor 2 d = 16 to 16/4 = 4. Hence, for d = 4, staggered fermions 
represent 4 flavors of mass-degenerate fermions. Staggered fermions have a single 1-component pair of 
Grassmann variables x x an d Xx per lattice point x. The corresponding lattice action for free staggered 
fermions takes the form 

S[x, X\ = a 4 IT (x x Vx,nXx+n ~ Xx+fJIx^Xx) +^Y1 m X x Xx, (47) 

x,fj, ^ a x 

where 

For m = the staggered fermion action has an exact U(l) e <S> U(1) Q symmetry 

X f x = exp(iip e )xx, X'x = X x exp(-iy? ), for (x 1 + x 2 + x 3 + x 4 )/a even, 
X' x = exp(iip a )xx, Xx = Xexp(-zy? e ), for (xi + x 2 + x 3 + x A )/a odd, (49) 

which is a subgroup of the SU(4)l ® SU(4:)r ®U(1)b chiral symmetry of the corresponding continuum 
theory. In the interacting theory the chiral and flavor symmetries besides U(l) e ® U(l) are explicitly 
broken by the staggered fermion action. These symmetries are recovered only in the continuum limit. 
A detailed introduction to staggered fermions and its properties can be found in textbooks like (23 |2H] 
and review articles like [40j which cover more advanced topics. 

Unlike Nature's QCD, staggered fermions have 4 flavors of the same mass. Since they have a rem- 
nant of chiral symmetry and are relatively easy to simulate numerically, staggered fermions provide a 
convenient framework for studies of chiral symmetry breaking at Nf = 4. Some numerical studies of 
staggered fermions try to mimic Nature's QCD by taking roots of the fermion determinant in order to 
relate 4 mass-degenerate staggered fermions to physical up, down, and strange quarks. These calcula- 
tions are not on completely solid grounds before one can show that locality (which may be violated at 
finite lattice spacing) is recovered in the continuum limit. 

3.5 Perfect Lattice Fermions 

In this subsection we relate the continuum theory of free fermions to a corresponding lattice theory by 
an exact renormalization group transformation. This is achieved by defining lattice fermion fields as 
block averages of continuum fields integrated over hypercubes IT3] . The resulting lattice theory 

is in all respects equivalent to the underlying continuum theory, i.e. it is completely free of lattice 
artifacts. For example, it has the same energy-momentum dispersion relation as the continuum theory. 
Even more important, it has an exact chiral symmetry (which may, however, be hidden). Lattice actions 
with these properties are known as perfect actions. 
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Let us derive a perfect fermion action by blocking from the continuum [Tz]. For this purpose we 
average the continuum fermion field ip(y) over hypercubes c x of size a d centered at the points x of a 
ci-dimensional Euclidean lattice 

^ X = -J d d y i>(y), ^ X = -J d d y ^(y), (50) 

a J c x CL J c x 

which in momentum space corresponds to 

*(p) = J2 *P(P + 2vr//a)n(p + 2-Kl/a), V(-p) = ^(~P - 2im/a)Tl{p + 2vm/a), (51) 

Note that the lattice fermion field is periodic over the Brillouin zone. The Fourier transform of the 
blocking kernel is given by 

= 2sin( P , fl /2) - (52) 

The lattice fermion propagator is related to the continuum propagator by 

(¥(-p)*(p)) = ^ (^(-p-2vr//a)^(p + 27r//a))n(j9 + 2vr//a) 2 
«ez d 

= E[^(^ + 2^/a)+m]- 1 n(p+27r//a) 2 . (53) 
zez d 

For m = the lattice propagator corresponds to a lattice action 

S[^,^] = a d Y / ^P^-y)^y, (54) 

with couplings p^x — y) calculable by a Fourier transformation. This lattice action is perfect by 
construction, i.e. its spectrum is identical with the one of the continuum theory. Hence, there should 
be no fermion doubling. On the other hand, the action is manifestly chirally invariant. This seems to 
contradict the Nielsen-Ninomiya theorem. However, the theorem is evaded because the action turns out 
to be nonlocal. Its couplings p M (x — y) do not decay exponentially at large distances. Instead for d>2 
they decay only power-like. As a consequence, in momentum space p M (p) is not regular (it actually 
has poles) and therefore the topological arguments behind the Nielsen-Ninomiya theorem do not apply. 
The nonlocality can be seen easily for d — 1. Then 71 = 1 and the sum in eq.([53|) can be performed 
analytically, resulting in a massless propagator that takes the form 

(*(-pMp)) = Y\i{p + 27r//a)]- 1 n(p + 2txI /a) 2 = £ cot(^). (55) 

This implies 

Pi(p) = -tan(^), (56) 
a 2 

which is singular at the edge of the Brillouin zone (p = ±7r/a). The corresponding coupling in coordinate 
space, 

Pl (x-y) = k-l) (x - y)/a , (57) 
a 

does not decay at all at large distances x — y and thus describes an extremely nonlocal action. For 
d > 2 the chirally symmetric perfect action remains nonlocal with a power-law decay of the couplings 
at large distances. 

Although the nonlocality of the perfect action arose naturally by blocking the theory from the 
continuum, from a practical point of view it is very inconvenient. For example, in a numerical simulation 
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it would be very demanding to include couplings to far away neighbors. It follows from the Nielsen- 
Ninomiya theorem that, in order to obtain a local perfect action, one must break chiral symmetry 
explicitly. This can be done by modifying the above way of blocking from the continuum which was 
chirally covariant. If one chooses to break chiral symmetry explicitly in the blocking procedure, the 
resulting perfect lattice action is not manifestly chirally invariant, but it is local. This can be achieved 
by constructing a perfect lattice action S[&, *$f] as 

exp(-S\9, *]) = J Vi)V^ exp{-^p J d d p ^(-p) [z 7/ ^ + m]V(p)} 

I I r _ 

x exp{ — -— - / d d p [¥(-p) - ]T ip(-p-27rn/a)Il(p + 2im/a)} 

x [tf (p) - V<P + 2vr//a)n(p + 2vrZ/a)]. (58) 
iez d 

The coefficient c is a source of explicit chiral symmetry breaking, which is injected into the theory via 
the renormalization group transformation that maps the continuum theory to the lattice theory. For 
c^Owe recover the chirally invariant but nonlocal perfect lattice action from before. In general one 
obtains 

(¥(-p)*(p)) = £ [Z7>„ + 2nya) + m]~ 1 n(p + 2nl/a) 2 + c, (59) 
zez d 

which corresponds to a local perfect action as long as c 7^ 0. 

Let us vary c in order to optimize the locality of the perfect action. For this purpose we again 
consider d — 1. Then the sum in eq. (|59|) can be performed analytically and the fermion propagator 
takes the form 

2 mo, . .pa., i 

-=— cothl ) — 1 coti — 

m m 2 a l y 2 ' v 2 

If we choose 

exp(ma) — 1 — ma 



(tf(-p)tf(p)) = f: - r ^[coth(^) -zcot(^)]" 1 + c. (60) 



(61) 



the propagator reduces to 

/ exp(ma) — 1 



(*(-p)*(p)> 



ma 



.1 . , . exp(ma) — 1 2 . 2 ,pa N 

sin (pa J H 1 — sin ( — 

a a a 2 



(62) 



This corresponds to the standard Wilson fermion action except that the mass m is now replaced by 
(exp(ma) — l)/a. Hence, for the above choice of c, in one dimension the perfect action is ultralocal, i.e. 
it has only nearest-neighbor interactions. In the massless limit m = the optimal choice for locality 
is c = a/2. When we go to more than one dimension the action remains local, but it is no longer 
ultralocal. 

Next we derive the energy-momentum dispersion relation of perfect lattice fermions. The fermion 
2-point function takes the form 

{¥(-/! 0)*(/r. .,-,,)) = ^ f W/a dp d (y(-p)V(p))exp(ip d x d ) 

[*7mG°m + 27r/ M /a) + m]~ 1 II(p + 2-kI/o) 2 + c > exp(ip d x d ) 
iez d J 

dp d ^2 ~ 

leZ d-i (p + 2-kI/ a) 2 +p d + m 2 
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^ ( 2sm{j> i a/2) \ 2 ( 2 sm(p d a/2) \ 2 

x 11 ^T~T exp{ip d x d ) + c S Xd>0 

f-Ji V Pi® + J \ Pda J 

= E C(p + 2irl/a)exp(-E(p + 27rl/a)x d ) + c5 Xdfl . (63) 
Fez** -1 

The sum over l d has been combined with the integral of p d over [— n/a, it /a] to an integral over the 
momentum space of the continuum theory. The sum over the spatial I G Z d_1 leads to infinitely many 
poles of the integrand, and hence to infinitely many states that contribute an exponential to the 2-point 
function. The energies of these states are given by the location of the poles, E(p + 2irl/a) = —ip d , with 

E(p + 2nt/a) 2 = -p\ = (p + 2ixl/a) 2 + m 2 . (64) 

Hence, the energy-momentum dispersion relation of perfect lattice fermions is exactly the same as in 
the continuum. In particular, as a result of exact blocking from the continuum, there are no lattice 
artifacts. Furthermore, the form of the renormalization group blocking transformation has no effect on 
the physical spectrum. In particular, the explicit chiral symmetry breaking term proportional to c only 
leads to a contact term c $ Xd ,o m the 2-point function. Hence, it has no effect on the spectrum which is 
extracted from the 2-point function at large Euclidean time separations x d . Remarkably, the spectrum 
of the lattice theory displays the consequences of Poincare invariance despite the fact that the lattice 
action only has the discrete lattice symmetries. 

Chiral symmetry is hidden in a similar way. Due to the explicit chiral symmetry breaking parameter 
c in the renormalization group blocking transformation, even for m = the perfect lattice action is not 
manifestly chirally invariant. Still, all physical consequences of chiral symmetry are correctly reproduced 
by the perfect action. As we will see later, this is due to the by now famous Ginsparg- Wilson relation 

{D-\ lb } = a lb . (65) 

Here D is the lattice Dirac operator and D~ l is the lattice fermion propagator. Indeed, using the 
optimal parameter c = a/2 for the perfect action one obtains 

{D-\ 75} = { E l^M + ^/a^Uip + 2nl/a) 2 + c, 75} = 2c 75 = «75- (66) 
iez d 

The Ginsparg- Wilson relation is the key to understanding chiral symmetry on the lattice [Hj. In the 
continuum, chiral symmetry implies {D~ l , 75} = 0. If one insists on this relation also on the lattice, i.e. 
if one insists on manifest chiral symmetry for a lattice action, the Nielsen-Ninomiya theorem implies 
fermion doubling (or, even worse, a violation of locality). The Ginsparg- Wilson relation {-D -1 , 75} = 075 
reduces to the relation {-D -1 , 75} = in the continuum limit a — > 0. Still, at finite lattice spacing a, the 
right-hand side of the Ginsparg- Wilson relation implies an explicit breaking of chiral symmetry. In the 
case of the perfect action the explicit breaking is due to the parameter c = a/2 in the renormalization 
group blocking transformation. This minimal explicit violation of chiral symmetry is sufficient to evade 
the Nielsen-Ninomiya theorem, and thus to avoid fermion doubling. Still, as we have seen explicitly for 
the perfect action, the physics (for example, the spectrum) remains the same as in the continuum. We 
will see later that the Ginsparg- Wilson relation leads to a natural definition of chiral symmetry on the 
lattice which reduces to the usual one in the continuum limit. 

4 Lattice QCD 

While many properties of lattice chiral symmetry can be studied in the free theory, one certainly also 
needs to understand it in the interacting theory. Hence, it is now time to endow the quarks with their 
nontrivial QCD dynamics by coupling them to the gluon field. However, before doing so, we first discuss 
lattice Yang-Mills theory without quarks. 
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4.1 Lattice Yang-Mills Theory 

Maintaining manifest gauge invariance is essential when gauge theories are regularized on the lattice. 
In the continuum, gauge transformations involve space-time derivatives of group- valued functions Q(x). 
On the lattice there are no infinitesimally close points, and continuum derivatives are usually simply 
replaced by finite differences. However, in order to maintain gauge invariance, one must proceed more 
carefully. Wegner and Wilson, as well as Smit, independently introduced the concept of a parallel 
transporter U x ^ G SU(N C ) connecting neighboring lattice points x and x + jl. The parallel transporter 
is related to an underlying continuum gauge field A^(x) = igA a ^(x)T a by 

U Xlll = Pexp f a dt A^x + fit), (67) 

<j 

where V denotes path-ordering. Under a non-Abelian gauge transformation the parallel transporter 
transforms as 

= n.tf.Xt*- (68) 

Wilson has constructed the Yang-Mills part of a simple lattice QCD action by multiplying parallel 
transporters around an elementary plaquette. The standard Wilson action is constructed as a sum over 
all plaquettes 

Sym[U] = -a 4 £ -l^Tr[U x ,,U x+fl , u Ul +0 ,,Ul u + U^U^M+fiAj- (69) 

x,fi,v 9 a 

It reduces to the continuum Yang-Mills action in the limit a — > 0. 

To fully define the path integral we must also consider the measure. The lattice functional integral 
is obtained as an integral over all configurations of parallel transporters U Xjft , i.e. 

Z = U[ du ^ exp(-S YM [U}). (70) 



X,fl 



One integrates independently over all link variables using the local Haar measure dU^ x for each parallel 
transporter. The Haar measure is a left- and right-invariant measure, i.e. 

dU f(ttU) = f dU f(Un) = f dU f(U), (71) 

SU(N C ) JSU(Nc) JSU(N C ) 

for any function f(U) and for any SU(N C ) matrix Q. It is convenient to normalize the measure such 
that 

/ dU=l. (72) 

JSU(N C ) 

For compact groups like SU(N C ) the integration is over a finite domain. This makes it unnecessary to 
fix the gauge in lattice QCD because the functional integral is finite even without gauge fixing. This is 
another important advantage of the formulation using parallel transporters. 

The Yang-Mills functional integral from above contains a single parameter — the bare gauge coupling 
g. When one wants to perform the continuum limit, one must search for values of g for which the 
correlation length of the lattice theory diverges in lattice units. In the language of statistical mechanics 
one is looking for a second order phase transition. Due to asymptotic freedom, in lattice QCD one 
expects a second order phase transition at g — > 0. To analyze the phase structure of a gauge theory one 
needs to study order parameters. A simple local order parameter like (U X)li ) is not useful. This follows 
from Elitzur's theorem jUj which states that gauge-variant observables simply vanish. A useful order 
parameter in a gauge theory must be gauge invariant and, in addition, nonlocal. In a pure gluon theory 
a good order parameter was suggested independently by Wegner and Wilson as 

W c = Tt I] U Xt(t . (73) 
(x,fi,)ec 



17 



For a rectangular curve C with side lengths R and T the Wilson loop describes the instantaneous 
creation and annihilation of a static quark-anti-quark pair at distance R which then exists for a time 
T. The Wilson loop is related to the static quark-anti-quark potential V(R) by 

lim (W c ) ~ exp(-V(R)T). (74) 

T — >oo 

In QCD we expect quarks and anti-quarks to be confined to one another by a potential rising linearly 
at large separations R, i.e. 

lim V(R) ~ <tR, (75) 

R^oo 

where a is the string tension. In a confinement phase the Wilson loop hence shows an area law 

lim (W c ) ~exp(-a J RT). (76) 

R,T — >oo 

Confinement is indeed verified very accurately in numerical simulations of lattice Yang-Mills theories. 
4.2 Standard Wilson Action for Lattice QCD 

We still need to couple the quarks to the gluons. First we do this by gauging SU(N C ) in the action of 
free Wilson fermions of eq.(j45j) 

SqcdI^, U] = a 4 £ o~(^*7AM**+/i ~ + a 4 E m^ x ^ x 

X,fJ, U X 

- fl4 E 4^^^+A^+i>,X" + U x , u Ux + ^Ul + ^UlJ- (77) 

X,fl,U 9 ^ 

In order to eliminate the doubler fermions we have introduced the Wilson term which breaks chiral 
symmetry explicitly. The lattice regularized functional integral takes the form 

Z= JJ f M x d9 x H f dU Xitl exp(-S QC D^, *,V))- (78) 



It depends on two parameters — the bare gauge coupling g and the bare quark mass m. Due to 
asymptotic freedom, in order to reach the continuum limit one must take g — > 0. When one puts m = 
for free Wilson fermions one reaches the chiral limit. In the interacting theory, on the other hand, this 
is no longer the case. In particular, since chiral symmetry is explicitly broken, the bare quark mass 
m must be fine-tuned in order to reach a massless limit. This fine-tuning is very unnatural from a 
theoretical point of view. In particular, following the discussion in the introduction, if one imagines 
Wilson's lattice QCD as an oversimplified model for the short- distance physics at the Planck scale, 
one could not understand at all why there are light fermions in Nature. The fine-tuning of m is also 
inconvenient from a practical point of view. For example, in a numerical simulation of Wilson's lattice 
QCD one must fine-tune m to many digits accuracy in order to make the pion massless. If one does 
this at relatively large g, i.e. before one reaches the continuum limit g — > 0, this massless "pion" is not 
even a proper Goldstone boson of a spontaneously broken chiral symmetry. For Wilson fermions an 
exact chiral symmetry that can break spontaneously does not exist at finite lattice spacing. It emerges 
only in the continuum limit after a delicate fine-tuning of m. 
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4.3 Ginsparg- Wilson Relation and Liischer's Lattice Chiral Symmetry 

In the discussion of the perfect free fermion action we have encountered the Ginsparg- Wilson relation 
eq. f|65|) . Following Liischer [21., we will now use this relation to construct a version of chiral symmetry 
that is natural for a lattice theory and reduces to the usual one in the continuum limit. For this purpose 
we consider a lattice fermion action 

% U] = a 4 ^D[U]9 = a 4 J2^ D PUA ( 79 ) 

x,y 

which is defined in terms of the lattice Dirac operator D[U]. This operator should be local (i.e. it should 
decay exponentially at large distances x — y) but not ultralocal. The lattice Dirac operator obeys the 
Ginsparg- Wilson relation if the corresponding fermion propagator D\U]~ l satisfies 

{D[U}-\ 75} = D[U]~ l l5 + j 5 D[U}- x = a 75 . (80) 

Alternatively, the Ginsparg- Wilson relation can be written as 

l5 D[U] + D[C7] 75 = aD[U] l5 D[U}. (81) 

It is nontrivial to construct lattice actions that obey the Ginsparg- Wilson relation. Until now we have 
seen that the perfect action for massless free fermions indeed satisfies this relation. In the next two 
subsections we will see that the same is true for perfect actions for massless QCD as well as for overlap 
fermions. For the moment we don't worry about the concrete form of Z?[J7], we just assume that it 
obeys eq.(jST|). 

Let us first consider an infinitesimal chiral rotation of the form familiar from the continuum 

= * + = (1 + ie a T a ^)m, 

= ^> + 5^> =^>{l + ie a T a ^). (82) 

Here T a (with a G {1, 2, iV| — 1}) are the generators of SU(Nf) and e a is a small parameter. In 
order to discuss flavor-singlet axial transformations with an infinitesimal parameter e° we also define 
T = 1. If the lattice action is local and has no fermion doubling, the Nielsen-Nino miya theorem 
implies that it cannot be invariant under the above chiral rotations. On the other hand, the lattice 
fermion measure is invariant under the full chiral symmetry U(Nf)i, <g> U(Nf)n. This is very different 
from massless QCD in the continuum. In the continuum the action is invariant under U (Nj)l® U (Nf)n 
chiral transformations, while the measure is invariant only under SU(Nf)i, ® SU(Nj)jt ® U(1)l=r- In 
particular, due to the anomaly the measure of the continuum theory is not invariant under flavor-singlet 
axial transformations, while the measure of the lattice theory is invariant. 

Next we consider Liischer's modification of the standard chiral transformation 

= + S^j = (1 + ie a T a l5 (l - ~D[U))\ % 

Tjf = *j7 + W = * ^1 + ie a T a (l - ^D[U])j 5 ^ . (83) 

Through D[U] Liischer's lattice version of a chiral transformation depends on the gluon field. Still, in 
the continuum limit a — > it reduces to the standard chiral symmetry of the continuum theory. It is 
remarkable that eq.([83p is a symmetry of any lattice action that obeys the Ginsparg- Wilson relation 
eq.([Blj). This follows from 

y'D[U]y' = * (l + ie a T a {l - ^D[U]) l5 ) D[U] (l + ie a T a l5 {l - ^D[U])) * 

= VD[U]V + V {ie a T a [ l5 D[U} + D[U]j 5 - aD[U]j 5 D[U}]) V + 0{e 2 ) 

= VD[U]ty + 0{e 2 ). (84) 
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Similarly, the variation of the lattice fermion measure takes the form 



Vm'W = ie a T a {l - ~£>[£/]h 5 ) (l - ^ a T a 75 (l - ^D[U]j) 

= VHW (1 - ie a Tr[T a l5 {2 - aD[U})]) + 0{e 2 ) 

= VVVV (l + ie°aTrfr 5 D[U]]) + 0(e 2 ). (85) 



Hence, while any Ginsparg- Wilson fermion action is invariant under Liischer's lattice chiral symmetry, 
the lattice fermion measure is not. Exactly as in the continuum, the fermionic measure of the lattice the- 
ory changes under flavor- singlet axial transformations, while it is invariant under SU(Nf)L®SU(Nf)p(3 
U{\)l=r- We will see later that the non-invariance of the fermionic measure under flavor-singlet axial 
transformations indeed gives rise to the correct axial anomaly. 



4.4 Classically Perfect Action for Lattice QCD 

Hasenfratz and Niedermayer have initiated and carried out an impressive program of explicitly con- 
structing nonperturbative lattice actions that are perfect at least at the classical level [J2|- The pure 
gauge part of a classically perfect action is the fixed point of a renormalization group blocking transfor- 
mation that maps lattice gauge fields [Zona fine lattice with lattice spacing a to lattice gauge fields U' 
on a coarser lattice with spacing 2a. Denoting the blocking kernel by Tym[U, U'], the classically perfect 
action obeys the condition 

S YM [U'} = nun (S YM [U] + T YM [U, U'}) , (86) 

which implicitly defines S[U]. This equation can be solved numerically by an iterative procedure on a 
multi-layer of finer and finer lattices. 

The fermionic part of a classically perfect massless fermion action results as a fixed point of a 
renormalization group transformation which maps a fermion field \I/ on the fine lattice to a fermion field 
\&' on the coarse lattice using a blocking kernel T F [U], i.e. 

= T F [U]V %, = ]T T F [U] X ,^ X . (87) 

X 

Here x is a point on the fine lattice and x' is a point on the coarse lattice. The fermionic part of the 
classically perfect action is quadratic in the fermion fields and can be written as 

S F $, % U] = a 4 UD P [U\^ = a 4 Y.^* D p\ U \^y ( 88 ) 

The corresponding fermion propagator Dp\U]~ 4 obeys the equation 

Dpp'Y 1 = T F [U]D P [U}- l T F [U} + c D P [U']-} y , =Y,T F [U} X , , x D P [U}- 4 y T F [U} y ^ + c S^. (89) 

x,y 

Here U is the gauge field on the fine lattice that minimizes the expression on the right-hand side of 
eq.(|BEJ) given a gauge field U' on the coarse lattice. The parameter c is analogous to the one introduced 
in the perfect action of the free theory. For an appropriate fermionic blocking kernel T F [U] and a 
choice of c that again optimizes locality one can show that the classically perfect action obeys the 
Ginsparg- Wilson relation. 

There are a number of important properties of the classically perfect fixed point action. For example, 

(i) The free gauge and fermion parts of the action reproduce the exact relativistic spectrum of the 
continuum. 
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(ii) It is possible to define a topological charge such that the action of any gauge field configuration 
with topological charge Q is larger than 87r 2 \Q\/ g 2 , a relation that is known from the continuum 
theory. 



(iii) There is an exact index theorem on the lattice that relates the fermionic zero modes of the fixed 
point Dirac operator and the topological charge of the lattice gauge field. 

(iv) Due to the Ginsparg- Wilson relation there is an exact chiral symmetry. 

For more details we refer the reader to the article by Hasenfratz JS] and references therein. 

Although classically perfect actions satisfy a number of attractive properties, they are difficult to 
implement in practice. A highly nontrivial and crucial step is to find practical (but still accurate) 
parameterizations of the perfect action. Although the initial work involved spin and gauge models 
in two dimensions, recently practical approximations of the fixed point actions for lattice QCD have 
emerged. In the pure gauge theory a useful parametrization was found in jTUj, while the fermionic 
problem was tackled in The resulting approximations of classically perfect actions yield physical 
results that are almost completely free of lattice artifacts |4"51 I44j . 

4.5 Domain Wall and Overlap Fermions 

In the early nineties Kaplan proposed a novel method to preserve chirality on the lattice The 
idea was to use the fact that chiral fermions become trapped on domain walls [3]. Kaplan used a 
Wilson-Dirac operator in five dimensions with a mass term that is a function of the fifth direction. In 
particular, the mass term changes sign creating a four-dimensional domain wall at the points where 
it vanishes. A four- dimensional chiral fermion is then trapped on the domain wall. In the meantime 
Narayanan and Neuberger were developing an idea of using an infinite number of regulator "flavor" 
fields to preserve chirality jSj. They realized that Kaplan's construction was equivalent to their idea 
since the fifth dimension is analogous to a flavor space. They used their interpretation and argued that 
the determinant of a chiral fermion in the background of a gauge field is equivalent to the overlap of two 
many-body fermionic ground states [8J. Initially, it seemed that the overlap was a reliable technique to 
regulate even chiral gauge theories on a lattice. Unfortunately, it was realized that the fermions with 
opposite chirality (which originate from an anti-wall) cannot be easily decoupled |1SJ US] • 

Although it was not clear whether the domain wall and overlap approach gave a completely satisfac- 
tory formulation of lattice chiral gauge theories involving no doublers of opposite chirality, something 
highly nontrivial had been achieved. It was possible to construct a lattice theory with a fermion with 
opposite chirality and use it effectively to preserve chiral symmetry in a vector-like gauge theory. An 
elegant way to use a five- dimensional fermion action to represent quarks in lattice QCD was first pro- 
posed by Shamir jl7j and elaborated further by Furman and Shamir in • This fermion is commonly 
referred to as a domain wall fermion and is used extensively in lattice simulations. Its action is con- 
structed on a five- dimensional space-time lattice with coordinates (x, 25), where x refers to the usual 
four dimensions and x^ G {05, 2as, L5} refers to the fifth direction of finite extent L5. Since the fifth 
direction is physically different from the other directions we have introduced a new lattice spacing 05 
in that direction. The domain wall fermion action is given by 



S F [V,V,U] =a 4 a 5 **,**D DW [V] 



x,x s ;y,y B 




(90) 



x,x 5 ,y,y5 



The domain wall Dirac operator is given by 
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{ {Pr S 2as , ya ~ <W)/a 5 - mP L <5 L5iy5 for x 5 = 1, 

D± [U]*b,V5 = \ ( P R S x 5 +a 5 ,y 5 + 1*1. <\r, a 5 ,y 5 - 4 5 ,S/5)/ a 5 for Gt 5 < S < L 5) (91) 

[ (Pl h 5 -a 5 ,y 5 - $L 5 ,y 5 )/a 5 - mP R 5 a5:y5 for x 5 = L 5 . 

Here Pr and P^ are the chiral projection operators defined in eq.(|18|). In the above action the parameter 
M is not the mass of the quark that is bound to the wall. By comparing with Wilson fermions one 
sees that the sign of the Wilson term has changed. In order to produce massless quarks one should 
set < Ma§ < 2 at tree level |5] and take L5 —>■ 00. There is a technical problem that needs to be 
taken into account. When L5 becomes infinite there are only Nj flavors of four-dimensional massless 
quarks bound to the wall, but there is an infinite number of modes at the cut-off. This may cause 
spurious effects at low energies. Hence, one needs to use bosonic (Pauli-Villars type) fields to cancel 
the contribution of these high-energy modes. For a detailed discussion of how this can be accomplished 
we refer the reader to |48j . 

There is a close connection between the domain wall approach and the overlap formula developed 
by Narayanan and Neuberger. Neuberger realized that it is possible to find an analytic formula for an 
effective Dirac operator that describes the massless chiral mode of the domain wall fermion. Using his 
insight on the overlap formula for vector-like gauge theories jlHl , he found a simple and elegant formula 
for the four- dimensional Dirac operator JS]) which is referred to as the overlap Dirac operator and 
which is given by 

D [U] = ^\l + l5 ^B=\ (92) 



2a 



H[U]< 



where H[U] = 7sP"[77] and P"[£7] is the operator we defined above in the context of the domain wall 
fermion. In order to obtain massless quarks one needs to set < Ma$ < 2 as before. In fact, it is 
possible to find an analytic formula for an effective Dirac operator that represents the chiral massless 
modes of the domain wall fermion even for finite L5 |49| IT^ E01 OH]- This operator takes the form 

Du [U] = i [1 + • //r:) ""' ( ' //rV " K , (93) 
J ' (l + H[U]) L */ a * + {l-H[U]) L */ a ** 



a 



where 



^=^M*m = 2r3m- (94) 



In the limit of L 5 — > 00 one obtains 

jE/l=£W [tfl = ^-[l + 75- 

>H[U\ 



lim D L5 [U] = D DWO [U] = ^[l + 75 Z I (95) 



l 5 ^oo -' 2a 

which reduces to the overlap Dirac operator Dq[U] when a 5 — > 0. 

Although today we know that the Ginsparg- Wilson relation leads to an exact chiral symmetry on 
the lattice, this connection was not appreciated until recently. After Ginsparg and Wilson discovered 
this interesting relation they found it difficult to explicitly construct a local operator that satisfies it. 
The relation was soon forgotten. After the discovery of perfect and overlap fermions the relation was 
rediscovered by Hasenfratz. It is straightforward to check that Do[U] and Ddwo[U] indeed satisfy the 
Ginsparg- Wilson relation. These new Dirac operators couple every pair of sites on the lattice. It is 
possible to show that if one wants to benefit from good chiral properties of Ginsparg- Wilson fermions, 
one has to give up the notion of ultralocal actions [22] • However, as has been shown in [SS] , close to the 
continuum limit the couplings in the overlap Dirac operator fall off exponentially with the distance. In 
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this sense these new Dirac operators are still local. Unfortunately, the closeness to the continuum limit 
is quite important to maintain both the chiral and locality properties of the Dirac operator. Recently, 
a physical picture based on the locality of zero modes of the Dirac operator was used to map out the 
regions in coupling constant space where Do[U] and Ddwo[U] lead to a good regularization of massless 
quarks 

There have been recent efforts to generalize the Ginsparg- Wilson relation and use this as a guide 
to construct new classes of Dirac operators [56J. Since Dirac operators which satisfy the Ginsparg- Wilson 
relation exactly are computationally very demanding, there have also been efforts to find perfect Dirac 
operators that satisfy the Ginsparg- Wilson fermions approximately Another approach has been to 
expand the most general lattice Dirac operator in a basis of simple operators. The coefficients of the 
expansion then are determined using the Ginsparg- Wilson relation 58J. This approach has been used 
to construct a practical operator for lattice simulations [59J. 

4.6 Quenched Approximation 

The chiral limit of lattice QCD is computationally very demanding because all algorithms contain a step 
in which the inverse of the Dirac operator in a fixed gauge field background needs to be computed. If the 
Dirac operator has small eigenvalues this step becomes very time consuming. Further, in the presence 
of small eigenvalues the number of inversions necessary before a statistically independent configuration 
is generated also increases leading to inefficiencies in the algorithms. 
The full partition function of lattice QCD is given by 

Z = J] f dV x dV x ]l I dU^ ex P (-S YM [U}) exp(-a d VD[U]y) 

x x,/i JSU(N C ) 

= II / dU *,» exv(-S YM [U\) detD[U}. (96) 

x,n JSU(Nc) 

Observables take the form 

(°) = \ II / dU *,» exp(-Sy M [£/]) detD[U] 0[U], (97) 

A X,fJ, JSU ( N c) 

where 0[U] is the observable in the background gauge field U. For example, when the observable is 
a fermion bilinear, 0[U] is constructed from a quark propagator. The difficulties in the chiral limit 
arise due to the determinant factors in the above relations. In order to avoid such difficulties, the 
determinant factor is sometimes dropped. This is an approximation in which one ignores the effects 
of virtual quark loops, which is commonly referred to as the quenched approximation. It should be 
pointed out that the quenched approximation is not a systematically controlled one. It is possible to 
take into account some effects of virtual quark loops by allowing the quarks in the loop to be heavier. 
This is the so-called partially quenched approximation. Many calculations today use either a quenched 
or partially quenched approximation to compute physical quantities for light quarks. 

Not surprisingly, in the chiral limit the quenched approximation introduces unphysical singularities. 
These singularities are again related to the small eigenvalues of the Dirac operator. However, now 
they enter observables through the quark propagator. Using quenched chiral perturbation theory, one 
can predict the form of these singularities in various observables close to the chiral limit |6TH I61j. For 
example, the chiral condensate is expected to diverge logarithmically with a dimensionless coefficient 5. 
Such spurious divergences are referred to as quenched chiral logarithms. Recently, several groups have 
evaluated 5 with results varying from about 0.03 to about 0.26. For further details we refer the reader 

to ng. 

The quenched approximation is particularly severe in the case of Wilson fermions close to the chiral 
limit. This is because, due to the absence of an exact chiral symmetry, the Wilson Dirac operator may 
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contain arbitrarily small eigenvalues when the quark mass is small even if it is non-zero. In the case 
of staggered fermions and Ginsparg- Wilson fermions one can show that the eigenvalues of the Dirac 
operator cannot be smaller than the quark mass, as expected in the continuum. This problem hinders 
simulations with dynamical light Wilson quarks. 

By adding a so-called twisted mass term to the two-flavor Wilson fermion action, it is possible to 
regulate the small eigenvalues of the Dirac operator This term does not change the continuum 

limit of the theory, since one can argue that — provided the continuum limit is taken carefully — the 
extra term is equivalent to a redefinition of the fermion field jHlj. This new approach to lattice QCD 
with an additional twisted mass term is becoming another practical approach to study the effects of 
light quarks in QCD, which has been tested in p3o] . 

5 Special Features of Ginsparg- Wilson Fermions 

We have already discussed in section [Ql how the Ginsparg- Wilson relation leads to a new realization of 
chiral symmetry in a finite lattice theory. This lattice chiral symmetry makes domain wall and overlap 
fermions, which were introduced in section 14. 5\ special in various ways. In this section we discuss some 
of the features of these Ginsparg- Wilson fermions that allow us to relate lattice quantities to continuum 
physical quantities more directly than it was possible before. 

5.1 Anomaly and Topological Charge on the Lattice 

It is well-known that the flavor- singlet chiral symmetry is anomalous in QCD. This can be shown in 
the continuum in essentially two ways: (i) By taking suitable care of the ultraviolet regulator in the 
diagrammatic approach one can show that the divergence of the flavor- singlet axial current is equal to 
the topological charge density of the gluon field; (ii) in the functional integral representation of QCD 
the fermionic integration measure is not invariant under flavor-singlet axial transformations when the 
Dirac operator in the background of a gluon field configuration has a non-zero index EH! ■ Before the 
discovery of the Ginsparg- Wilson relation and its consequences, in a lattice formulation the anomaly 
could only be understood after a complex calculation [HE]- Since the fermionic measure was always 
invariant under chiral transformations, either the anomalous symmetry was explicitly broken by the 
lattice action (as in the case of Wilson or staggered fermions) or the fermion doubling would cancel the 
anomaly completely (as for naive fermions). The only viable approach was to calculate the appropriate 
diagrams in the lattice theory, which is usually quite tedious, and then to take the continuum limit. 

The discovery of the exact chiral symmetry of Ginsparg- Wilson fermions allows one to derive the 
anomaly in a straightforward fashion as explained in [21] . Here we review the essential steps of that 
proof. Consider the expectation value of an operator Ofi&, \&] in the background of a gauge field U 
given by 



where D[U] XtV is a lattice Dirac operator that obeys the Ginsparg- Wilson relation. For simplicity we 
have suppressed color, flavor, and Dirac indices. Using the variation of the fermionic measure eq.([85|) 
under Luscher's infinitesimal flavor-singlet axial transformations as well as the invariance of the action 
one can show that 






(99) 



a(z - D[U]) l5 (z - D[U}) = z{2 - az) l5 - (1 - az)[(z - D[U]) l5 + 75^ - D[U\)\, 



(100) 
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where z is a complex number not contained in the spectrum of D[U]. After multiplying both sides of 
this equation from the right with (z — D\U])~ l and taking the trace one gets 

- aTr[ l5 D[U]] = z(2 - az)Tr[ l5 (z - Dp})- 1 }. (101) 

We can now divide both sides by the factor z(2 — az) and integrate over a small circle centered at the 
origin that does not encircle any spectral value of D[U] other than 0. Since 

Po = f^~ (z-DlU])- 1 (102) 

projects on the subspace of zero modes of D[U], we then get 

aTr[-f 5 D[U]} = 2 (n_ - n+) = 2 index(L> [[/]). (103) 

where n± represent the number of zero modes of D[U} which are also eigenstates of 75 with eigenvalues 
±1. This yields the anomalous Ward-identity 

(SO) = 2 mdex{D[U])(0), (104) 

which is familiar from the continuum but which is now valid in a completely regularized finite lattice 
theory. 

Often a lattice Dirac operator D[U] is 75-Hermitean, i.e. it obeys the relation = 75D75, in 
addition to satisfying the Ginsparg- Wilson relation. In that case one can simplify the above discussion 
which relates aTr^Dft/]] to the index of -D[i/]. This was done in |oTH ITU] where the lattice index 
theorem was first derived. One can show that the eigenstates u\ of aD[U] are such that the eigenvalues 
A fall on a circle given by A = 1 — e ia . There are three types of eigenstates: (i) Those with A = 
such that 75MA = ±«a- Let n± represent the number of such eigenstates with eigenvalue ±1; (ii) those 
with A = 2 such that 75WA = ±«a- Let n' ± represent the number of such eigenstates with eigenvalue 
±1; (iii) those with complex (non-real) A, in which case 75WA — Since Tr(7 5 ) = 0, we must have 
n_ — n + = n' + — n'_. Using these relations it is easy to derive eq. ([103)1 . 

The above calculation shows that there is indeed an anomaly in a theory with Ginsparg- Wilson 
fermions. However, it is also important to show that the anomaly reproduces the well-known results 
in the continuum limit. This has been shown perturbatively for a class of overlap Dirac operators in 
[7T] and nonperturbatively for the overlap Dirac operator of eq.(!92|) in [721 EH EH EH! • More recently, 
the nonperturbative arguments were extended to the overlap Dirac operator that is constructed using 
a perfect fermion action . Thus, from a variety of studies we can conclude with confidence that the 
above calculation does indeed reproduce the correct anomaly in the continuum limit. 

In the continuum, the index of the Dirac operator is related to the topological charge of the back- 
ground gauge field configuration. On the lattice there is no unique way to define the topological charge. 
A naive discretization of the topological charge does not even yield an integer on the lattice. A geomet- 
ric construction of the topological charge is possible which then leads to integer values [77\. The index 
of a Ginsparg- Wilson Dirac operator is yet another definition of the topological charge. In this case we 
define the topological charge as 

Q=±aTt[ 7s D[V]]. (105) 

This definition is more meaningful since it is the same charge that enters the anomalous Ward identity 
in eq. (|104|) . As we will see in the next section, the Witten-Veneziano mass formula can be derived on 
the lattice using this definition of the topological charge. Further discussions of the anomaly and the 
topological charge can be found in a recent review [70] and the references therein. 
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5.2 Witten-Veneziano Mass Formula on the Lattice 



The mass of the ?/-meson is closely connected to the anomalous flavor-singlet chiral symmetry of QCD. 
In the absence of the anomaly one can argue that the flavor- singlet chiral symmetry would be sponta- 
neously broken leading to a massless Goldstone boson for massless quarks. The absence of a low-mass 
flavor- singlet pseudo-scalar particle in Nature with the correct flavor quantum numbers is usually at- 
tributed to the anomaly. To compute this mass from first principles is one of the outstanding challenges 
in QCD. Due to the nonperturbative nature of the physics involved, the only viable approach to this 
problem from first principles is using lattice QCD. 

Interestingly, the physics of the 7/-meson mass can be understood more easily in certain limits of 
QCD. For example, in the limit N c — > oo, with g 2 N c and Nf held fixed [7S], or by assuming that the 
anomalous flavor-singlet axial Ward identities retain their validity order by order in an expansion in 
u = Nf/N c around u = [75] . In both cases one can derive a leading-order Witten-Veneziano relation 
for the mass of the r/-meson which is given by 

2iV r 

m l = ~pk I ^ x (i( x M q ))ymi (106) 
where F n is the pion decay constant and q(x) is the topological charge density 

q{x) = -^W-Trt^)^)]- (107) 

Here the trace is over color indices. The subscript YM in eq. (jlU6|) indicates that the gg-correlation 
function is to be computed in the pure Yang-Mills theory, i.e. in the absence of quarks. 

Although the Witten-Veneziano relation is valid only in certain limits of QCD, one can use it 
to estimate m v i in a lattice QCD calculation. Unfortunately, eq. (jl06p is rather formal and cannot 
be translated to lattice QCD without addressing a number of subtleties (HOj 1 - Two main problems 
need to be solved in order to make eq. (jl()6)) rigorous and of practical use. One has to: (i) find a 
properly normalized lattice definition of the topological charge density q(x); (ii) subtract from q{x)q(0) 
appropriate contact terms, so as to define it properly in the continuum limit. As we have discussed 
earlier, different lattice definitions of q(x) are possible. Furthermore, the second problem can be quite 
subtle since the subtraction may involve contact terms of the form c 5(x) which contribute a finite 
constant to the r/-meson mass. Some of these subtleties were first pointed out in [8*2*] . The reader is 
also referred to a recent discussion in [£%] . 

One of the special features of Ginsparg- Wilson fermions is that we can use eq. (jl05|) to define the 
topological charge density q x (which obeys a 4 J2 X 1x = Q) as 

q x = ^Tr[j 5 D[U] a>x ). (108) 

In this expression the trace is only over the color and Dirac indices, but not over the space-time index 
x. As explained in [HI], this definition allows one to derive 



2 2iV> 4 



J2(^Qx)ym, (109) 



in the appropriate limit in lattice QCD. Thus one obtains a relation for the 7/'-meson mass exactly as in 
the continuum. Thanks to the chiral symmetry properties of Ginsparg- Wilson fermions, no additional 
subtractions are necessary. 



1 For a recent review we refer the reader to [ST] and the references therein. Here we review only some basic points. 
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5.3 Renormalization of Operators 

Computing the matrix elements of fermionic currents in QCD is necessary for determining a variety 
of physical quantities such as hadronic decay constants, electromagnetic and weak form factors, quark 
masses, etc. A knowledge of the lattice renormalization constants of these operators is necessary to 
relate the matrix elements computed using lattice simulations to the corresponding ones defined in 
continuum renormalization schemes, like the MS scheme, used in experimental data analysis. Recently, 
several groups have contributed to this subject. We refer the reader to |B3| and references therein for 
a recent review of renormalization in lattice field theory. The existence of exact symmetries simplifies 
the computations of these renormalization constants enormously. Hasenfratz has discussed operator 
renormalization in the context of perfect actions [H3] . In this section, using a few examples, we illustrate 
how one can obtain constraints among various renormalization constants. A more general method to 
renormalize massless Ginsparg- Wilson fermions in lattice gauge theories has been worked out in [86J. 
For an extensive discussion of perturbative renormalization, especially in the context of Ginsparg- Wilson 
fermions we refer the reader to the recent review jSTj and references therein. 

Consider the renormalization constants for local bilinear quark operators of the form 

o a { x ) =${x)r a i>{x), (no) 

where T a denotes generic Dirac matrices, i.e. 1, 75, 7 M , 7^75, o^. Specific bilinear operators, denoted 
according to their Lorentz group transformations, are 

S{x) = ${x)i/>(x), P{x) = ^(x) l5 i/;(x), Vfa) = ^(x) 7m T>(x), A%{x) = ^{x) W T a ^{x). (Ill) 

The lattice renormalization constants can be obtained by the equation 

r?W) = ^M^Mr^a), (112) 

where Tq(j>, a) and Tq S (p,h) are the two-quark one-particle irreducible vertex functions with an in- 
sertion of O calculated respectively on the lattice and in the continuum using the MS renormalization 
scheme, and Z^(afi) is the wave function renormalization constant. In the continuum, chiral symmetry 
imposes constraints among the various renormalization constants. For example, it is well-known that 
Z s = Zp and Zy = Such relations reduce the amount of work in computing the renormalization 
constants. On the lattice one often has to deal with mixing of operators with lower-dimensional oper- 
ators leading to additional subtractions that may diverge with an inverse power of the lattice spacing. 
Fortunately, with Ginsparg- Wilson fermions the renormalization of operators is also quite analogous to 
the continuum. For example, one can argue that the relations Zs = Zp and Zy = Za remain valid [88J. 

In order to see the usefulness of the Ginsparg- Wilson relation, let us define a new set of fermion 
bilinear operators 

o' a = *r a (l - ^D[U\) ¥, (113) 

where D[U] is the Ginsparg- Wilson Dirac operator appearing in the fermionic action of the theory. The 
addition of the term depending on D[U] is harmless since its contributions vanish in the continuum 
limit. This leads to the fact that the renormalization constants Zo> for the new operators are the same 
as Zq. In fact, these new "primed" operators are more natural in a Ginsparg- Wilson theory at finite 
lattice spacings. As pointed out in jHS] and (HO], for finite lattice spacing the spontaneous breaking of 
non-singlet chiral symmetry is related to the vacuum expectation value of the operator 5", i.e. it occurs 
if (S') 0. Further, by using the infinitesimal flavor-singlet chiral transformations given in eq.([83|) one 
can show that 

8S'(x) = 2s°P'(x), 5P'(x) = 2e°S'(x). (114) 
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Under the non-singlet chiral transformations the vector and axial currents V^ a and A'^ transform as 

5v ,a = iftcjA* 5A >« = ifteJv*. (115) 

Using the corresponding Ward identities one can show that 

Zgi = Zp/, Zyi = Zj\i. (116) 

Thus, we find that Z$ = Zp and Zy = Za- 

There is also a simplification in determining the renormalization of the quark mass using Ginsparg- 
Wilson fermions when the bare fermion mass is introduced in the action by writing the massive Dirac 
operator as 

D m {U] = D[U] +m(l- \d[U]\ ■ (117) 

The renormalization of the fermion mass is then related to that of the operator 5" denoted by Z$- In 
particular, it is easy to show that 

Z m = Z s \ (118) 

The renormalized quark mass is then obtained as m ren = Z m m. 

Using these ideas, one-loop perturbative computations of these renormalization constants were per- 
formed with overlap fermions in [SH] , while a nonperturbative approach has been used in jUTj to deter- 
mine the quark condensate and quark masses. Nonperturbative renormalization has also been studied 
extensively with domain wall fermions, the details of which can be found in [52*] . 



5.4 Lattice Simulations with Ginsparg- Wilson Fermions 

Over the past few years, numerical simulations with overlap fermions, domain wall fermions, and ap- 
proximately perfect fermions have begun to appear. Computationally, these fermions are about fifty 
to a hundred times more expensive than conventional fermions. Thus, the most reliable results from 
lattice simulations using these new fermions are still obtained in the quenched approximation where 
the fermion determinant is ignored. In this section we review some of the recent work. 

Simulations with domain wall fermions are in much better shape thanks to both a larger number of 
researchers working on this subject and faster computers being used in the computations. A detailed 
study of the chiral properties of domain wall quarks can be found in [53] and [HI]- Recently, an extensive 
calculation of quenched QCD with light quarks using overlap fermions was performed in [95J. Chiral 
properties of overlap fermions were studied earlier in |DTH 157] . The first results using approximately 
classically perfect fermions, which obey the Ginsparg- Wilson relation to a good accuracy, have also 
been published in |55| . 

Calculations using domain wall fermions typically use lattices of size 24 3 x 40 to 32 3 x 60 at a lattice 
cut-off which varies from 1/a = 2 GeV to 3 GeV. On the other hand, simulations with overlap fermions 
are currently done on lattices of sizes from 12 3 x 24, 16 3 x 32 to 20 4 at lattice cut-offs of 1/a = 1.33 
GeV to 2 GeV. Results using approximately perfect fermions have been obtained on 8 3 x 24, 12 3 x 24, 
and 16 3 x 32 lattices at lattice cut-offs 1/a = 1.3 GeV, 2 GeV, and 2.5 GeV. 

A variety of quantities can be calculated in lattice simulations. This includes the hadron spectrum, 
the topological susceptibility, the quark condensate, the pion decay constant, the average of the up and 
down quark masses m = (m u + md)/2, the strange quark mass m s , weak matrix elements like the kaon 
Bjc parameter, properties of the nucleon, such as its axial charge qa and the parameters that arise in 
quenched and unquenched chiral perturbation theory. A comparison of the results obtained in different 
studies indicates that there are systematic errors in many of these quantities that are presently of the 
order of 10 — 20%. For example, the topological susceptibility defined as 

* - M (119) 
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is found to vary from about (0.176GeV) 4 in j^S], (0.196GeV) 4 in [H, to (0.213GeV) 4 in |96j. Finite 
temperature effects were recently studied in |lUUj and it was found that Xt varies from about (0.191GeV) 4 
at T = 0.88T C to (O.lOOGeV) 4 at T = 1.31T C where T c is the deconfinement temperature. Currently, 
several groups are performing a careful analysis of the systematic errors due to finite volume effects, 
finite lattice spacing errors, as well as the effects of dynamical quark loops. 

Another important parameter in QCD is the chiral condensate. However, unlike Xt this quantity 
needs a proper definition in a renormalization scheme |lUlj . For Ginsparg- Wilson fermions a nonper- 
turbative renormalization group invariant definition of the condensate was introduced in jUTj. Using 
this definition the condensate was found to be (0.243(10)GeV) 3 in |lU2j using overlap fermions and 
(0.235(ll)GeV) 3 in |lU3j using the perfect action approach. When these are converted to the MS 
scheme at a scale of 2GeV one obtains (0.266(15)GeV) 3 , which is consistent with (0.250(3)GeV) 3 ob- 
tained in |95j . These values are in agreement with the value (0.256(8)GeV) 3 obtained in j^Hl using 
domain wall fermions. The errors in the chiral condensate seem to be at the 4 — 8% level. 

One of the most impressive achievements in lattice QCD is the determination of the hadron spectrum. 
It is indeed exciting that one can compute the masses of a variety of hadrons from first principles and 
compare them with experiments. In the quenched approximation this spectrum shows about 10% 
deviations from the experimental results [104J. These deviations appear to reduce when one introduces 
dynamical quarks |lU5j . Since most calculations are not close to the continuum limit, to some extent the 
deviations depend on the quantity used to determine the lattice spacing. Recently, results using chirally 
improved fermions in the quenched limit have begun to emerge. Results using fixed point fermions 
suggest that deviations from the experimental results are strongly correlated with their experimentally 
observed widths |98| . In other words, unstable particles are more severely affected by quenching than 
stable ones. Using domain wall fermions, the quenched spectrum has been investigated in 106J. 

Several groups have recently obtained results on kaon physics using overlap and domain wall 
fermions. The main goal in these projects is to extract the kaon Bk parameter, which plays a central 
role in understanding CP violation in the kaon system. The value of Bk depends on the renormalization 
scale and scheme. When evaluated at fi — 2 GeV in the MS scheme the value appears to be around 
0.57 with an error of about 5 — 10%. We refer the reader to the work in |107t 11081 1109] IllOj for further 
details. Quark masses have also been studied in [^3 llllj using overlap fermions and in |109t I112j using 
domain wall fermions. The current estimates are m = 3.5 — 4 MeV and m s = 99 — 133 MeV. 

Recently, domain wall quarks have been used to compute the nucleon axial charge ■ One such 
calculation finds = 1.21(5) in the chiral limit. This should be compared with the experimental 
result gA = 1.267(3). We refer the reader to |113j for further details and for a comparison with earlier 
work. There are also interesting results on the excited states of the nucleon from different types of 
chiral fermion actions |114[ 11151 1116] . Finally, we like to mention that efforts to compute the parton 
distribution functions from lattice QCD are well-developed (see |117j for a recent review). Recently, 
a method to measure even generalized parton distribution functions has been proposed |118j and has 
been applied using domain wall fermions |119j . Such calculations have also been performed using Wilson 
fermions jUHl QUI H22J • 

There are many other studies we have not touched on in this article. All these studies indicate that, 
with sufficient computing power, at least in the quenched approximation, computations using Ginsparg- 
Wilson fermions are feasible. Dynamical fermion calculations are also being envisioned. However, in 
the absence of an algorithmic breakthrough for dynamical fermions, it is likely that most computations 
will be limited to quark masses that are still too large for chiral perturbation theory to work reliably. 
This may keep us from studying the full dynamics of chiral symmetry from first principles for some 
time. 
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6 Low-Energy Effective Theories 



Due to spontaneous chiral symmetry breaking, at low energies the QCD dynamics is dominated by 
Goldstone bosons [123J — the pions for Nf = 2. It is possible to use a low-energy effective description 
that only involves the Goldstone boson fields |124| I125| I126j . Chiral perturbation theory provides a 
systematic low-energy expansion that predicts the pion dynamics based on symmetry principles and 
a few low-energy parameters (like the pion decay constant F n and the chiral order parameter (ipip)) 
whose values can be determined either from experiments or from lattice QCD calculations |127j . The 
low-energy interactions of pions and nucleons can be understood in baryon chiral perturbation theory 
1281 11291 11301 H31| I132j . Effective theories that are relevant to more than one nucleon can be useful 
even in the absence of pions |133| H34j. For a recent review of this subject we refer the reader to [135J. 

Usually, low-energy effective theories are treated perturbatively directly in the continuum. In order 
to investigate nonperturbative effects that may arise within the effective theory, it is interesting to also 
regularize it on the lattice 1.3(3 H37| I138j . There are efforts underway to understand nuclear matter on 
the lattice |139| 11401 H41| starting from some effective theory. Interestingly, the nonlinear realization of 
chiral symmetry on the lattice does not lead to the same subtleties (like the fermion doubling problem) 
that one faces within the microscopic QCD theory |142j . 



6.1 Effective Theory for Goldstone Bosons 

The Goldstone bosons are the lightest particles in QCD. Therefore they dominate the dynamics of the 
strong interactions at low energies. It is possible to use a low-energy effective description that only 
involves the Goldstone boson fields. This is not only true for QCD but also for any other system with a 
continuous global symmetry G which is spontaneously broken to a subgroup H. The Goldstone bosons 
are described by fields in the coset space G/H in which points are identified if they are connected by 
symmetry transformations of the unbroken subgroup H. In QCD we have G/H = SU(Nf). Hence the 
Goldstone boson fields are represented by special unitary matrices U(x) G SU(Nf). Under global chiral 
rotations they transform as 

U(x)' = LU(x)R ] . (120) 

Goldstone bosons interact weakly at low energies. Their effective Lagrangian is constructed as a deriva- 
tive expansion. The leading term of the pion effective action takes the form 



S[U] = J d 4 x 



4 Zl\ f 



:i2D 



The first term on the right-hand side is chirally invariant. Its prefactor is the pion decay constant F n 
which determines the strength of the interaction between the Goldstone bosons. The second term is 
the chiral symmetry breaking mass term which contains the quark mass matrix. Under chiral transfor- 
mations this term transforms as 

Tr{MU' j + U'PV) = Tr(MRtftf + LUWM ] ). (122) 

If all quark masses are equal, i.e. if M. = ml, the Lagrangian is invariant against SU(Nf)p flavor 
rotations for which R = L. For a general diagonal mass matrix the flavor symmetry is reduced to 

ll/VU!;, 

The constants F n and (ipip) determine the low-energy dynamics at leading order and enter the effec- 
tive theory as a priori unknown parameters. These parameters can be determined from experiments or 
through nonperturbative lattice QCD calculations. Up to these two low-energy constants the Goldstone 
boson dynamics is completely determined by chiral symmetry. At higher energies additional terms arise 
in the effective theory. Again, they are restricted by chiral symmetry and they contain new low-energy 
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parameters — the Gasser-Leutwyler coefficients. Chiral perturbation theory (with mass-degenerate 
quarks) is a systematic low-energy expansion around the classical vacuum configuration U(x) = 1. One 
writes 

U(x) = exp(2m a (x)T a /F 7r ), a E {1, 2, Nj-1}, (123) 

where T a are the generators of SU(Nj), and one then expands in powers of n a (x). In this way one can 
derive, for example, the Gellmann-Oakes-Renner relation 



6.2 Effective Theory for Nucleons and Pions 

Chiral perturbation theory can be extended to sectors with non-zero baryon number. Nucleons enter 
the low-energy effective theory in the form of a Dirac spinor field ip(x) and ip(x) that transforms as an 
SU(2)j isospin doublet. Global chiral rotations L <E> R € SU(2)l <g> SU(2)r can be realized nonlinearly 
on this field using the transformations 

rj}{x)' = V(x)ip(x), i)(x)' = ^(x)V(xf. (125) 

The field V(x) depends on L and R as well as on the field U(x) and can be written as 

V(x) = R(R^LU(x)) 1/2 (U(x) 1/ y = L(L^RU(xf) 1/2 U(x) 1/2 . (126) 

For transformations in the unbroken isospin vector subgroup SU(2)j = SU(2)l=r the field V(x) reduces 
to the global flavor transformation V(x) = L = R. The local transformation V(x) is a nonlinear 
representation of chiral symmetry which has the form of a local SU(2) transformation, despite the fact 
that it represents just a global SU{2) L <g> SU{2) R symmetry. 

In order to construct a chirally invariant (i.e. SU (2) "gauge" invariant) action one needs an SU(2) 
flavor "gauge" field. For this purpose one constructs a field u(x) e SU(2) from the pion field U(x) as 

u(x) = U(x) 1/2 . (127) 

The matrix u(x) is located in the middle of the shortest geodesic connecting U(x) with the unit-matrix 
11 in the group manifold of SU(2). Under chiral rotations the field u(x) transforms as 

u(x)' = Lu(x)V(xf = V(x)u(x)Rl (128) 

The anti-Hermitean composite field 

v n( x ) = 7i[ u ( x V^u(x) + u(x)d (l u(x)\ (129) 



transforms as a "gauge" field 

M x )' = \[V(x)u(xfL^d^Lu(x)V(xf) + V(x)u(x)rfd^Ru(xfV(xf)] = V(x)(v^(x) + d^)V(xf. 

(130) 



A Hermitean composite field is given by 



a n( x ) = -\u{xfd^u(x) - u(x)dnu(xy], (131) 



which transforms as 



a^x)' = \[V{x)u{x)^L^d^Lu{x)V{x)^) - V(x)u{x)R^d fjL {Ru{xfV{xf )} = V(x)a^{x)V{xf . (132) 
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The leading terms in the Euclidean action of a low-energy effective theory for nucleons and pions 
take the form 

S[ip, ip,U] = J d 4 x {M^i) + ipJuidft + v^ip + igA^Pl^a^ 

+ S-Trld^d^U] - ^-Ti[MU ] + M ] U}}. (133) 

Here M is the nucleon mass generated by spontaneous chiral symmetry breaking and gA is the coupling 
to the isovector axial current. It is remarkable that — thanks to the nonlinear realization of chiral 
symmetry — the fermion mass term is chirally invariant. This makes sense because the mass M arises 
from the spontaneous breakdown of chiral symmetry even in the chiral limit. It is remarkable that 
fermions with a nonlinearly realized chiral symmetry do not contribute to anomalies \V28\ I143j . In the 
low-energy effective theory anomalies enter through the Wess-Zumino-Witten term |144|. 1145] . 



6.3 Effective Theory for Constituent Quarks 

The chiral quark model of Georgi and Manohar |146j is formulated in terms of gluons, pions, and 
constituent quarks ip{x) and ip(x) which transform in the fundamental representations of SU(N C ) and 
SU(Nf). In particular, under the nonlinearly realized SU(Nf)L <8> SU(Nf)n chiral symmetry the con- 
stituent quark field transforms as 

iP( x y = V(x)ip(x), ${x)' = ^V(x)\ (134) 
while under an SU(N C ) color gauge transformation 

ip(x)' = Q(x)ip(x), ${x)' = ^n(x)l (135) 
The Euclidean action of the chiral quark model is given by 

S[ip, ip, U,A] = J d A x {Mipip + ilrfpidp + + A^ip + igA^l^a^ 

+ ^Trldprfd^U] - ±(fal>)Tr[Mrf + M^U] - ^TrfF^F^]}. (136) 

In this case, M is the constituent quark mass which is generated by spontaneous chiral symmetry 
breaking. It should be noted that in the effective theory the strong gauge coupling g is weaker than in 
QCD itself. However, due to the coupling to the gluon field the constituent quarks are still confined. 

The chiral quark model is based on the assumption that the energy scale for chiral symmetry breaking 
is larger than the one for confinement. An effective description in terms of constituent quarks which 
receive their mass from chiral symmetry breaking before they get confined by residual low-energy gluons 
should then make sense. In fact, the phenomenological success of the nonrelativistic quark model may 
suggest that this picture is indeed correct. In their work, Georgi and Manohar provided a framework 
that puts the nonrelativistic quark model on a solid field theoretical basis. 

A potential problem of the chiral quark model is related to the confinement scale. The value of g 
may be so small that there are unacceptably low-lying glueball states. For the same reason, there might 
be a low-temperature deconfinement phase transition in the gluon sector significantly below the finite 
temperature chiral phase transition. These potential problems are impossible to address quantitatively 
in the continuum formulation of the chiral quark model because they involve the nonperturbative 
dynamics of confinement. In order to be able to address these issues, it is useful to formulate the 
constituent quark model on the lattice. 
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6.4 Lattice Formulation of Nonlinearly Realized Chiral Symmetry 

Let us construct theories with a nonlinearly realized chiral symmetry on the lattice. The Goldstone 
boson field U x G SU(Nf) naturally lives on the sites a; of a four-dimensional hypercubic lattice and it 
transforms as 

U' x = LU X R\ (137) 

under global chiral rotations. As in the continuum, the field u x G SU(Nf) is constructed as u x = U^ 2 
which transforms as 

Lu x V} = V x u x Rl (138) 



We now proceed to the construction of the lattice analog V X:fl G SU(Nf) of the continuum flavor 
"gauge" field v^(x) which is a flavor parallel transporter along a lattice link in the group SU(Nf). In 
analogy to the continuum expression eq. ()129|) for v^(x) we construct 



V*,p = Y a ^ Ux+ > 1 + u * u Uii} ■ ( 139 ) 

In the continuum limit V XtlM = exp[av^(x + /t/2)]. However, at finite lattice spacing V xli is in general 
not an element of SU(Nf), just a complex Nf x Nf matrix in the group GL(Nf). One can project a 
group- valued parallel transporter V x>fi G SU(Nf) out of V Xifl by performing a GL{Nf) / SU(Nf) coset de- 
composition |142j . By construction, under the nonlinearly realized chiral symmetry V x ^ then transforms 
as a parallel transporter, i.e. 

= v*v^vU- ( 14 °) 

Next we construct a lattice version of the continuum field a^(x) defined in 1)131)) . For this purpose, 
we first construct 

While in the continuum a M (x)' = V{x)a IJi {x)V{x)\ the lattice field A x>fl transforms as 

Also, in contrast to the continuum field a M (x), the lattice field A x ^ is in general neither traceless nor 
Hermitean. It is therefore more natural to introduce the field 

4^ = \\K,Vl, + V^AlJ - ^Ltt[A x ,,V^ + V Xifl AtJl, (143) 

which, by construction, is traceless and Hermitean and which transforms as 

A L J = V X A^V X \ (144) 
with the matrix V x located at the site x on the left end of the link (x, u). Similarly, we define the object 

At = \K,A X ^ + A^V X J - ^Tr[Vl,A x>fl + A^V X J1, (145) 

which transforms as 

K! = ( 146 ) 

with the matrix Vx+f, located at the site x + fi on the right end of the link (x, /i) . It should be noted 
that A x „ and A x „ are not independent but are related by parallel transport, i.e. 

Ax,n = Vl^V^. (147) 
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6.5 Constituent Quarks on the Lattice 



Using the lattice construction of a nonlinearly realized chiral symmetry presented in the previous sub- 
section, it is now straightforward to put, for example, Georgi and Manohar's chiral quark model on the 
lattice. The resulting fermion action takes the form 

+ <r £ ^(2M7.,.vi/,. - m7, i ;,,„/ -,.,/!/.,. . - - w r . ; ,i ,.). (148) 

x,n /a 

Here U XyfJ- G SU(N C ) denotes the standard Wilson color parallel transporters living on the lattice links. 
It is interesting to ask if the gluon dynamics of QCD can be modeled successfully in the chiral quark 
model. Since the chiral quark model does not represent a systematic low-energy expansion of QCD, 
one should not expect to obtain quantitative results directly relevant to QCD. However, even if only 
qualitative insight into the success of the nonrelativistic quark model can be gained, this would be quite 
interesting. 

The Wilson term in the above action removes the doubler fermions. As we discussed before, in 
standard lattice QCD this term breaks chiral symmetry explicitly. Remarkably, when chiral symmetry 
is nonlinearly realized, not only the fermion mass term (proportional to M) but also the Wilson term 
is chirally invariant. The only source of explicit chiral symmetry breaking is the current quark mass 
matrix hA. It is interesting to ask how the Nielsen- Ninomiya theorem jSS] has been avoided. Clearly, 
the action of eq. (|148j) is local. The corresponding Dirac operator is given by 

D[V,A, U] x>y = 
+ 
+ 



E 7^(24,, - VrJ-,^,-,, - \'r „.„)• (149) 



The Nielsen- Ninomiya theorem assumes that the Dirac operator anti-commutes with 75. This is not 
the case when chiral symmetry is nonlinearly realized, and hence one of the basic assumptions of the 
Nielsen-Ninomiya theorem is not satisfied. Interestingly, Ginsparg- Wilson fermions evade the Nielsen- 
Ninomiya theorem by violating the same assumption. Of course, the nonlinear realization of chiral 
symmetry requires an explicit pion field which is not present in the fundamental QCD Lagrangian. 

Remarkably, not only the lattice fermion action but also the lattice fermion measure is gauge invari- 
ant. At first sight this seems to be a severe problem because gauge invariance of both the fermion action 
and the fermion measure implies that lattice fermions with nonlinearly realized chiral symmetry do not 
contribute to anomalies. One might even suspect that the doubler fermions have not been properly re- 
moved and thus have canceled the anomalies of the physical fermions. Fortunately, this is not the case. 
Indeed fermions with a nonlinearly realized chiral symmetry do not contribute to anomalies. Instead, as 
mentioned in the continuum discussion, the anomalies are contained in the Wess-Zumino-Witten term 
|144[ 1145] which must be added to the lattice action explicitly. 
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7 Conclusions 



In this review we have discussed some basic issues of chiral symmetry on the lattice. The recent lattice 
developments have put chiral symmetry on a solid theoretical basis at a nonperturbative level. As 
discussed in the introduction, lattice QCD can now explain nonperturbatively why nucleons can exist 
naturally, i.e. without fine-tuning, far below the Planck scale, provided that space-time has additional 
hidden dimensions. This in turn explains why gravity is so weak, a nontrivial result one obtains from 
lattice QCD without doing any numerical work. Even chiral gauge theories like the standard model 
have now been constructed rigorously beyond perturbation theory. This is a very substantial step 
forward in the theoretical formulation of the basic laws of Nature. The development of Ginsparg- Wilson 
lattice fermions is also beginning to revolutionize practical lattice QCD simulations. In particular, if 
new algorithmic developments go hand in hand with the recent theoretical insights, Ginsparg- Wilson 
fermions may lead to substantial progress towards an accurate numerical solution of QCD. As usual, 
many new questions arise based on the new insights. For example, supersymmetry still waits to be put 
on rigorously solid grounds beyond perturbation theory. Also many practical numerical calculations 
with Ginsparg- Wilson fermions still need to be done. We hope that we have provided a certain basis 
for newcomers to enter this very active field of current research. 



Acknowledgements 

Over the years we have discussed chiral symmetry both in the continuum and on the lattice with many 
colleagues including O. Bar, T. Bhattacharya, W. Bietenholz, R. Brower, G. Colangelo, N. Christ, 
E. Farhi, J. Gasser, M. Gockeler, J. Goldstone, M. Golterman, P. Hasenfratz, R. Jackiw, R. Jaffe, 
J. Jersak, D. Kaplan, A. Kronfeld, M. Laursen, H. Leutwyler, M. Liischer, A. Manohar, T. Mehen, 
J. Negele, R. Narayanan, H. Neuberger, F. Niedermayer, M. Pepe, G. Schierholz, Y. Shamir, S. Sharpe, 
J. Smit, R. Springer, F. Steffen, and F. Wilczek. We gratefully acknowledge the insights they generously 
shared with us. This work was supported in part by funds provided by the U.S. Department of En- 
ergy (D.O.E.) under cooperative research agreement DE-FG02-96ER40945 and by the Schweizerischer 
Nationalfond (SNF). 

References 

[1] F. Wilczek, arXiv:hep-ph/0201222| 

[2] K. G. Wilson, Phys. Rev. D 10 (1974) 2445 

[3] K. G. Wilson, New Phenomena In Subnuclear Physics, ed. A. Zichichi, (Plenum Press, New York), 
Part A (1977) 69 

[4] C. G. Callan and J. A. Harvey, Nucl. Phys. B 250 (1985) 427 

[5] D. B. Kaplan, Phys. Lett. B 288 (1992) 342 

[6] R. Narayanan and H. Neuberger, Phys. Lett. B 302 (1993) 62 

[7] R. Narayanan and H. Neuberger, Phys. Rev. Lett. 71 (1993) 3251 

[8] R. Narayanan and H. Neuberger, Nucl. Phys. B 412 (1994) 574 

[9] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B 414 (1994) 785 



35 



[10] T. DeGrand, A. Hasenfratz, P. Hasenfratz, and F. Niedermayer, Nucl. Phys. B 454 (1995) 587 

[11] U.-J. Wiese, Phys. Lett. B 315 (1993) 417 

[12] W. Bietenholz and U.-J. Wiese, Nucl. Phys. B 464 (1996) 319 

[13] W. Bietenholz, R. Brower, S. Chandrasekharan, and U.-J. Wiese, Nucl. Phys. B 495 (1997) 285 
[14] P. Hasenfratz, Nucl. Phys. B 525 (1998) 401 

[15] P. Hasenfratz, S. Hauswirth, K. Holland, T. Jorg, F. Niedermayer, and U. Wenger, Int. J. Mod. 
Phys. C12 (2001) 691 

[16] P. Hasenfratz, Nucl. Phys. B (Proc. Suppl.) 63 (1998) 53 

[17] P. H. Ginsparg and K. G. Wilson, Phys. Rev. D 25 (1982) 2649 

[18] H. Neuberger, Phys. Lett. B 417 (1998) 141 

[19] H. Neuberger, Phys. Rev. D 57 (1998) 5417 

[20] H. Neuberger, Phys. Lett. B 427 (1998) 353 

[21] M. Liischer, Phys. Lett. B 428 (1998) 342 

[22] M. Liischer, Nucl. Phys. B 549 (1999) 295 ; Nucl. Phys. B 568 (2000) 162 
[23] A. Hasenfratz, Nucl. Phys. B (Proc. Suppl.) 119 (2003) 131 
[24] M. Liischer, JEEP 305 (2003) 052 
[25] D. B. Kaplan, ar Xiv:hep-lat/0 309099 

[26] M. Creutz, Quarks, Gluons and Lattices, Cambridge Monographs on Mathematical Physics, Cam- 
bridge University Press, (1983) 

[27] I. Montvay and G. Minister, Quantum Fields on a Lattice, Cambridge Monographs On Mathe- 
matical Physics, Cambridge University Press, (1994) 

[28] H. J. Rothe, Lattice Gauge Theories: An Introduction, World Scientific Lecture Notes in Physics, 
59 (1997) 

[29] R. Gupta, Introduction to lattice QCD, published in Les Houches, "Probing the standard model 
of particle interactions", Part 2 (1997) 83 

[30] M. Liischer, Advanced lattice QCD, published in Les Houches, "Probing the standard model of 
particle interactions", Part 2 (1997) 229 

[31] M. Creutz, Rev. Mod. Phys. 73 (2001) 119 

[32] H. Neuberger, Ann. Rev. Nucl. Part. Set. 51 (2001) 23 

[33] M. Salmhofer and E. Seiler, Lett. Math. Phys. 21, 13 (1991). 

[34] M. Salmhofer and E. Seiler, Commun. Math. Phys. 139, 395 (1991) [Erratum-ibid. 146, 637 
(1992)]. 

36 



[35] H. B. Nielsen and M. Ninomiya, Phys. Lett. B 105 (1981) 219 ; Nucl. Phys. B 185 (1981) 20 
[36] H. R. Quinn and M. Weinstein, Phys. Rev. Lett. 57 (1986) 2617 
[37] L. Susskind, Phys. Rev. D 16 (1977) 3031 

[38] H. S. Sharatchandra, H. J. Thun, and P. Weisz, Nucl. Phys. B 192 (1981) 205 

[39] H. Kluberg-Stern, A. Morel, O. Napoly, and B. Petersson, Nucl. Phys. B 220 (1983) 447 

[40] G. W. Kilcup and S. R. Sharpe, Nucl. Phys. B 283 (1987) 493 

[41] S. Elitzur, Phys. Rev. D 12 (1975) 3978 

[42] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B 414 (1994) 785 

[43] T. DeGrand, A. Hasenfratz, P. Hasenfratz, and F. Niedermayer, Nucl. Phys. B 454 (1995) 615 ; 
Phys. Lett. B 365 (1996) 233 

[44] P. Hasenfratz, S. Hauswirth, T. Jorg, F. Niedermayer, and K. Holland, Nucl. Phys. B 643 (2002) 
280 

[45] M. F. L. Golterman, K. Jansen, D. N. Petcher, and J. C. Vink, Phys. Rev. D 49 (1994) 1606 

[46] M. F. L. Golterman and Y. Shamir, Phys. Rev. D 51 (1995) 3026 

[47] Y. Shamir, Nucl. Phys. B 406 (1993) 90 

[48] V. Furman and Y. Shamir, Nucl. Phys. B 439 (1995) 54 

[49] R. Narayanan and H. Neuberger, Nucl. Phys. B 443 (1995) 305 



[50] Y. Kikukawa and T. Noguchi, |arXiv:hep-lat/ 9902022 
[51] L. Giusti, Nucl. Phys. B (Proc. Suppl.) 119 (2003) 149 
[52] I. Horvath, Phys. Rev. Lett. 81 (1998) 4063 

[53] P. Hernandez, K. Jansen, and M. Luscher, Nucl. Phys. B 552 (1999) 363 

[54] M. Golterman and Y. Shamir, Phys. Rev. D 68 (2003) 074501 

[55] K. Fujikawa, Nucl. Phys. B 589 (2000) 487 

[56] K. Fujikawa and M. Ishibashi, Nucl. Phys. B 605 (2001) 365 

[57] W. Bietenholz, Eur. Phys. J. C 6 (1999) 537 

[58] C. Gattringer, Phys. Rev. D 63 (2001) 114501 

[59] C. Gattringer, I. Hip and C. B. Lang, Nucl. Phys. B 597 (2001) 451 
[60] C. W. Bernard and M. F. L. Golterman, Phys. Rev. D 46 (1992) 853 
[61] C. W. Bernard and M. F. L. Golterman, Phys. Rev. D 49 (1994) 486 
[62] H. Wittig, Nucl. Phys. B (Proc. Suppl.) 119 (2003) 59 

37 



[63] W. A. Bardeen, A. Duncan, E. Eichten and H. Thacker, Phys. Rev. D 59 (1999) 014507 
[64] R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz [Alpha collaboration], JEEP 0108 (2001) 058 
[65] K. Jansen, A. Shindler, C. Urbach and I. Wetzorke [XLF Collaboration], |arXiv:hep-lat/0312013| 
[66] K. Fujikawa, Phys. Rev. Lett. 42 (1979) 1195 

[67] K. Fujikawa, Phys. Rev. D 21 (1980) 2848 [Erratum: Phys. Rev. D 22 (1980) 1499 ] 

[68] L. H. Karsten and J. Smit, Nucl. Phys. B 144 (1978) 536 

[69] P. Hasenfratz, V. Laliena, and F. Niedermayer, Phys. Lett. B 427 (1998) 125 

[70] F. Niedermayer, Nucl. Phys. B (Proc. Suppl.) 73 (1999) 105 

[71] T. Reisz and H. J. Rothe, Phys. Lett. B 455 (1999) 256 

[72] Y. Kikukawa and A. Yamada, Phys. Lett. B 448 (1999) 265 

[73] K. Fujikawa, Nucl. Phys. B 546 (1999) 480 

[74] H. Suzuki, Prog. Theor. Phys. 102 (1999) 141 

[75] D. H. Adams, Ann. Phys. (N. Y.) 296 (2002) 131 



[76] D. H. Adams and W. Bietenholz, |arXiv:hep- lat/0307022 
[77] M. Liischer, Commun. Math. Phys. 85 (1982) 39 
[78] E. Witten, Nucl. Phys. B 156 (1979) 269 
[79] G. Veneziano, Nucl. Phys. B 159 (1979) 213 

[80] B. Alles, M. D'Elia, and A. Di Giacomo, Nucl. Phys. B 494 (1997) 281 

51] L. Giusti, G. C. Rossi, M. Testa, and G. Veneziano, Nucl. Phys. B 628 (2002) 234 

[82] E. Seiler and I.O. Stamatescu, "Some remarks on the Witten- Veneziano formula for the eta' 
mass", preprint MPI-PAE/PTh 10/87 

[83] E. Seiler, Phys. Lett. B 525 (2002) 355 

[84] S. Sint, Nucl. Phys. B (Proc. Suppl.) 94 (01) 79 

[85] P. Hasenfratz, Nucl. Phys. B 525 (1998) 401 

[86] T. Reisz and H. J. Rothe, Nucl. Phys. B 575 (2000) 255 

17] S. Capitani, Phys. Rep. 382 (2003) 113 

C. Alexandrou, E. Follana, H. Panagopoulos, and E. Vicari, Nucl. Phys. B 580 (2000) 394 

[89] S. Chandrasekharan, Phys. Rev. D 60 (1999) 074503 

[90] Y. Kikukawa and A. Yamada, Nucl. Phys. B 547 (1999) 413 

[91] P. Hernandez, K. Jansen, L. Lellouch, and H. Wittig, JEEP 0107 (2001) 018 

38 



[92 
[93 
[94 
[95 
[96 
[97 

[98 
[99 

100 
101 
102 
103 

104 
105 

106 
107 
108 
109 
110 
111 

112 
113 

114 
115 



T. Blum et ai, Phys. Rev. D 66 (2002) 014504 



T. Blum et al, arXiv:hep- lat/0007038 

A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D 63 (2001) 114504 
T. W. Chiu and T. H. Hsieh, Nucl. Phys. B 673 (2003) 217 

T. DeGrand and U. M. Heller [MILC collaboration], Phys. Rev. D 65 (2002) 114501 

S. J. Dong, T. Draper, I. Horvath, F. X. Lee, K. F. Liu, and J. B. Zhang, Phys. Rev. D 65 (2002) 
054507 



C. Gattringer et al. [BGR Collaboration], |arXiv:hep-lat/0307013 



P. Hasenfratz, S. Hauswirth, K. Holland, T. Jorg and F. Niedermayer, Nucl. Phys. B (Proc. 
Suppl.) 106 (2002) 751 |arXiv:hep-lat/0109007| . 

C. Gattringer, R. Hoffmann and S. Schaefer, Phys. Lett. B 535 (2002) 358 
P. Hernandez, K. Jansen and L. Lellouch, Phys. Lett. B 469 (1999) 198 

P. Hernandez, K. Jansen, L. Lellouch and H. Wittig, Nucl. Phys. B (Proc. Suppl.) 106 (2002) 766 

P. Hasenfratz, S. Hauswirth, T. Jorg, F. Niedermayer and K. Holland, Nucl. Phys. B 643 (2002) 
280 

S. Aoki et al. [CP-PACS Collaboration], Phys. Rev. Lett. 84 (2000) 238 

A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. Lett. 85 (2000) 4674 ; [Erratum Phys. 
Rev. Lett. 90 (2003) 029902 ] 



Y. Aoki et al, |arXiv:hepdat/0211023 



T. DeGrand [MILC Collaboration], |arXiv:hep-lat/0309026| 

N. Garron, L. Giusti, C. Hobbling, L. Lellouch, and C. Rebbi, arXiv:hep-ph/0306295 
A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D 64 (2001) 114506 
T. Blum et al. [RBC Collaboration] , |arXiv: hep-lat /0 1 1 0075] 

L. Giusti, C. Hdlbling, and C. Rebbi, Phys. Rev. D 64 (2001) 114508 [Erratum: Phys. Rev. D 65 
(2002) 079903 ] 

T. Blum, A. Soni, and M. Wingate, Phys. Rev. D 60 (1999) 114507 

S. Sasaki, K. Orginos, S. Ohta, and T. Blum [RBC-KEK Collaboration], Phys. Rev. D 68 (2003) 
054509 

S. Sasaki, T. Blum, and S. Ohta, Phys. Rev. D 65 (2002) 074503 

D. Brommel, P. Crompton, C. Gattringer, L. Y. Glozman, C. B. Lang, S. Schaefer, and A. Schafer 
[BGR Collaboration] , |arXiv:hep-"ph/0307073l 



[116] S. J. Dong, T. Draper, I. Horvath, F. X. Lee, K. F. Liu, N. Mathur, and J. B. Zhang, 
|arXiv:hep-ph/0 306199 



39 



[117] J. W. Negele, Nucl. Phys. A 711 (2002) 281 

[118] P. Hagler, J. Negele, D. B. Renner, W. Schroers, T. Lippert, and K. Schilling [LHPC collabora- 
tion], Phys. Rev. D 68 (2003) 034505 

[119] W. Schroers et al. [LHPC collaboration] , |arXiv:hep-lat/0309065l 

[120] M. Gockeler et al, Phys. Rev. D53 (1996) 2317 

[121] M. G6ckeler et al, Phys. Rev. D63 (2001) 074506. 

[122] M. Gockeler et al, Phys. Rev. Lett. 92 (2004) 042002. 

[123] J. Goldstone, Nuovo Cimento A 19 (1961) 154 

[124] S. Coleman, J. Wess, and B. Zumino, Phys. Rev. 177 (1969) 2239 

[125] C. G. Callan, S. Coleman, J. Wess, and B. Zumino, Phys. Rev. 177 (1969) 2247 

[126] S. Weinberg, Physica 96A (1979) 327 

[127] J. Gasser and H. Leutwyler, Ann. Phys. (N.Y.) 158 (1984) 142 ; Nucl. Phys. B 250 (1985) 465 

[128] H. Georgi, Weak Interactions and Modern Particle Theory, Benjamin-Cummings Publishing Com- 
pany, (1984) 

[129] J. Gasser, M. E. Sainio, and A. Svarc, Nucl. Phys. B 307 (1988) 779 
[130] E. Jenkins and A. Manohar, Phys. Lett. B 255 (1991) 558 

[131] V. Bernard, N. Kaiser, J. Kambor, and U.-G. Meissner, Nucl. Phys. B 388 (1992) 315 
[132] T. Becher and H. Leutwyler, Eur. Phys. J. C 9 (1999) 643 

[133] D. B. Kaplan, M. J. Savage, and M. B. Wise, Phys. Lett. B 424 (1998) 390 ; Nucl. Phys. B 534 
(1998) 329 

[134] P. F. Bedaque, H.-W. Hammer, and U. van Kolck, Phys. Rev. C 58 (1998) 641 ; Nucl. Phys. A 
676 (2000) 357 

[135] P. F. Bedaque and U. van Kolck, Ann. Rev. Nucl. Part. Set. 52 (2002) 339 
[136] S. Myint and C. Rebbi, Nucl. Phys. B 421 (1994) 241 
[137] I. A. Shushpanov and A. V. Smilga, Phys. Rev. D 59 (1999) 054013 
[138] R. Lewis and P.-P. A. Ouimet, Phys. Rev. D 64 (2001) 034005 

[139] H. M. Miiller, S. E. Koonin, R. Seki, and U. van Kolck, Phys. Rev. C 61 (2000) 044320 
[140] J.-W. Chen and D. B. Kaplan, |aTXiv:hep-lat/0308016| 
[141] D. J. Lee and I. C. F. Ipsen, | |arXiv:nucl-th/0308052| 



[142] S. Chandrasekharan, M. Pepe, F. D. Steffen, and U.-J. Wiese, larXiv:hep-lat7 0306020 
[143] A. Manohar and G. W. Moore, Nucl. Phys. B 243 (1984) 55 

40 



[144] J. Wess and B. Zumino, Phys. Lett. B 37 (1971) 95 

[145] E. Witten, Nucl. Phys. B 223 (1983) 422 ; Nucl. Phys. B 223 (1983) 433 

[146] H. Georgi and A. Manohar, Nucl. Phys. B 234 (1984) 189 



41 



